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Abstract. 

We present a comprehensive analysis of the stellar population properties (age, metallicity and the alpha-element enhancement 
[E/Fe]) and morphologies of red-sequence galaxies in 24 clusters and groups from z ~ 0.75 to z ~ 0.45. The dataset, consisting 
of 215 spectra drawn from the ESO Distant Cluster Survey, constitutes the largest spectroscopic sample at these redshifts 
for which such an analysis has been conducted. Analysis reveals that the evolution of the stellar population properties of 
red-sequence galaxies depend on their mass: while the properties of most massive are well described by passive evolution 
and high-redshift formation, the less massive galaxies require a more extended star formation history. We show that these 
scenarios reproduce the index-cr relations as well as the galaxy colours. The two main results of this work are (1) the evolution 
of the line-strength indices for the red-sequence galaxies can be reproduced if 40% of the galaxies with cr < 175 kms"' 
entered the red-sequence between z = 0.75 to z = 0.45, in agreement with the fraction derived in studies of the luminosity 
functions, and (2) the percentage the red-sequence galaxies exhibiting early-type morphologies (E and SO) decreases by 20% 
from z = 0.75 to z = 0.45. This can be understood if the red-sequence gets populated at later times with disc galaxies whose 
star formation has been quenched. We conclude that the processes quenching star formation do not necessarily produce a 
simultaneous morphological transformation of the galaxies entering the red-sequence. 



Key words. Galaxies: high-redshift - Galaxies: stellar content - Galaxies: evolution - Galaxies: elliptical and lenticular - 
Galaxies: abundances 



1. Introduction 



Early - type galaxies are strongly clustered (e.g . . iLoveday et al 



1995 1: iHermit et al.lll996l: rWillmer et al.lll998l: IShepherd et al 



200 1*). making galaxy clusters an ideal place for their study. 
Stellar population studies have shown that both, local and high- 
redshift early-type galaxies follow tight correlations between 
the colors and line-strength features and other properties of the 
galax i es (mo stly related with their ma sses) (e .g. ^ower et al 



2005 



19921; iKuntsc hner 2000; Trager et alj l2000al: iBemardi et al. 
Sanchez-Blazquez et al.]2006bl) The mere existency of 



these correlations, as well as their evolution with redshift, 
seem to be consistent with a very early (z>2) and coordi- 
nated formation of their stars ('Visvanathan & Sandage"l977'; 
Bower et al. 1992; E ll is et al .199 7; Koda ma & Ari moto 1998; 
Stanford et al.1 1 19981; iKelson et all [200 ll; iBender et al.l 1 19981 
Ziegler et al. booih . However, the evolution of the cluster red- 



sequence luminosity function with redshift challenges this view 
(De Lucia et al. 2004; De Lucia et al. 2007; Kodama et al. 2004 
Rudnick et al. 2008, in preparation; although see AndreonI 



( 120081) for contradictory results), and so does the mo rphologi- 
cal evolution since z - 0.4 (e.g., iDressler et al.ll 19971) . and the 
evolution of the blue/star-forming fraction (Butcher & Oemler 
1984; Poggianti et al. 2006) in clusters. 

A deeper understanding of the nature of the evolution of 
the cluster red-sequence requires to go beyond colors (af- 
fected by the age-metallicity degeneracy) and derive the stel- 
lar populations parameters (age and chemical abundances) 
with time, using stellar population models. However, breaking 
the age-metallicity degeneracy requires a combination of in- 
dices with different sensitivities to both parameters (see, e.g., 
Rabin 1982). Q Unfortunately, instrumental limitations have 
long prevented accurate measurements of absorption line in- 
dices at high redshift. Recent observational advances, on the 
ground and in space, now allow such measurements in inter- 
medi ate and high redshi ft galaxies , both in clusters and in the 
field dZiegler et al. 200ll: lBarr et aL(2005l:(Sch"iavon et al.l2006t 



J0rgensen et al.a20 05). The analysis of these data have revealed 



that, indeed, the age-metallicity degeneracy is confusing the in- 
terpretation of the scaling relations: large spreads in galaxy lu- 
minosity weighted ages and metallicities at high redshift have 
been found in datasets showing very tight Faber-Jackson, Mgb- 
cr and Fundamental plane (FP) relations, challenging the clas- 
sical interpretation of the tightness of the scaling relations. 

High redshift spectroscopic samples are, however, still re- 
stricted to a maximum of ~30 galaxies. Moreover, they of- 



ten target a single cluster ('Kelson et al."2001^, 'j0rg ensen et al 



[200 5; Tran et al. 2007; Kelson et al. 2006). This could be bias- 
ing the results, as studies at low redshift have shown that the 
star formation histories of early-type galaxies might depend on 
cluster properties (e.g., compare Kuntschner & Davies 1998; 



* Based on observations obtained at the ESO Very Large Telescope 
(VLT) as part of the Large Programme 166.A-0162 (the ESO Distant 
Cluster Survey). 

' Combination of optical and near-IR colours have also been probed 
useful to break the age-metallicity degeneracy (e.g., Peletier, Valentijn 
(& Jameson 1990; Mac Arthur et al. 2004; James et al. 2005, among 
others), although dust reddening is still a problem. 



Caldwell et al. 20031: Sanchez-Blazquez et al. 2003 : J0rgensen 



19991: iNelan et al.ll2005l;lTrager et al.ll2008l) . Therefore, a study 
based on a large sample of galaxies in clusters covering a large 
range in cluster masses is necessary to obtain a complete pic- 
ture of galaxy evolution. 

The present work is based on 24 clusters with redshifts be- 
tween 0.39 and 0.8 from the ESO Distant Cluster Survey (here- 
after, EDisCS). The major novelty of the present work is that 
we do not only consider very massive structures; our sample 
spans a large range in cluster velocity dispersions. Therefore, 
we minimize possible biases due to the relationship between 
galaxy properties and environment. 

It is now common to study red-sequence galaxies 
rather than morphologically classified early-type galaxies, 
as colours are easier to measure in large datasets than 
morphology, and they appear to be r tiore correlated with 



environment dK auffma nn et alJ l2004t iBlanton et al.l 12005 



Martinez & Murielii2006l) . Nonetheless, observations have also 
revealed the existence o f an intrinsic spread in morphol- 
ogy at given colo ur (e.g.. IConselicell2006t iBaldry et alj|2004 : 
Cross et al ]|2004'). This stress the need for investigating the re- 
lation between morphology and colour, and their evolution with 
redshift and, therefore, this is the strategy we adopt. 

Our sample of red-sequence galaxies is much larger than 
those used in previous efforts. It encompasses 337 galaxies, 
distributed in 24 clusters and groups, for which we study stel- 
lar populations and morphologies. Our sample of red-sequence 
galaxies is not restricted to the most massive, but span a wide 
range of internal velocity dipsersion (100-350 km/s), compara- 
ble to the samples analysed at low-redshift. 

Throughout the paper, we adopt a concordance cosmology 
with Qm = 0.3, Qa = 0.7, Ho=70 km s"' Mpc"'. All magni- 
tudes are quoted in the Vega system. 

2. The sample 

The ESO Distant Cluster Survey (EDisCS) is a photomet- 
ric and spectroscopic survey of galaxies in 20 fields con- 
taining clusters with redshifts between 0.39 and 0.96. These 
fields w ere selected from the Las Campanas Distant Cluster 
Survey dGonzalez et al.ll200lh . specifically from the 30 high- 
est surface brightness candidate s. A full descripti on of the 
sample selection can be found in IWhite et alj (l2005h . EDisCS 
includes structures with velocity dispersions from ~150 to 
~1 100 km s"', i.e. from small groups to clusters. 



Deep optical photometry with VLT/FORS2 ( IWhite et al 



2005D, near-IR photometry with SOFI on the ESO/NTT 
(Aragon-Salamanca et al., in preparation), and deep-multi-sht 
spectroscopy with VLT/FORS2 were acquired for each field. 
The same high-efficiency grism was used in all observing runs 
(grism 600RIH-19, - 6780A ). The wavelength range varies 
with the field and the x-loc ation of the slit on the mask (see 
Milvang- Jensen et alj|2008l for details), but it was chosen to 
cover, at least, a rest-frame wavelength range from 3670 to 
4150 A (in order to include [Oil] and H6 lines) for the assumed 
cluster redshift. 

The spectral resolution is ~6 A (FWHM), corresponding to 
rest-frame 3.3 A at z = 0.8 and 4.3 A at z = 0.4. Typically, 
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four- and two-hour exposures were obtained for the high-z 
and mid-z samples, respectively. The spectroscopic selection, 
obser vations, data reduction and cata logs are presen ted in de- 
tail in lHanidaveT al. (2004) and Mil vang- Jensen et a l. (2008). 
In brief, standard reduction procedures (bias subtraction, flat 
fielding, cosmic -ray removal, geometrical distortion correc- 
tions, wavelength calibration, sky subtraction and flux cali- 
bration) were performed with IRAF. Particular attention was 
paid to the sky-subtraction that was performe d before any in- 



terpo lation or rebinning of the data (see Milva ng- Jensen et al 



20081 for details). This dataset has also been complemented 
wifli 80 orbits of H ST/ACS imaging in F814W of the high- 
est redshift clusters jPesai et al.ll2007n . Ha narrow-band imag- 
ing dPinn et al.ll2005b and X MM-Newton/EPIC X-ray obser 
vations ( iJohnson et al ] l2006l) . 

26 structures ( groups or clusters) were identified in the 



20 EDisCS fields (IHaUiday et al J 12004 iMilvang- Jensen et al 



12008 ). Two of these structures were not considered in this pa- 
per The first one, clll03. 7-1245, because its spectroscopic 
redshift (0.96) is too far away from the redshift targeted by 
the photo-z-based selection (0.70). This can introduce observa- 
tional biases. The second one, cll238.5-l 144, because it could 
only be observed for 20 minutes and the S/N of the spectra is 
not sufficient for our analysis. 

To define the red-sequence, we use total-I magnitudes, es- 
timated by adding, to the Kron magnitudes, a correction appro- 
priate for a point source measured within an aperture equal to 
the galaxy's Kron aperture. These corrections are obtained em- 
pirically from unsaturated and isolated stars on each convolved 
I-band image (see White et al. 2005 for details). 

A galaxy is considered a member of a given cluster 
(or group) if its redshift falls within Sxcrdustei of the clus- 
ter redshift, where crHmter i s the clust er velocity dispersion 



presen ted in iHalliday et alJ (120041) and iMilvang- Jensen et al 



(l2008h . Similarly to all other EDisCS analyses, galaxy groups 
are identified by crdustei < 400 km/s. We define red-sequence 
galaxies as those with secure redshifts, and colours between 
+0.3 mag from the best linear fit (with slope fixed to -0.09) 
to the colour-magnitude relation (V-I vs. I) of the objects 
without emission lines. This definition coincides with that of 



White et all (120051) and De Lucia et al. (2004, 2007). The width 
of the region corresponds to ~3 times the RMS -dispersion of 
the red-sequence colours of our two most populated clusters, 
C11216.8-1201 and cll232.5-1250. Table □ lists the structures 
selected for the present work and the total number of galaxies 
on the red-sequence. Figure [1] displays the V - / vs / colour- 
magnitude diagrams of the 24 EDisCS structures studied in the 
present paper. 

To study the evolution of the galaxy population properties 
with redshift, we divide our galaxies in 3 different redshift in- 
tervals : (1) 0.39 < z < 0.5; (2) 0.5 < z < 0.65; and (3) 
0.65 < z < 0.81. For convenience, we refer to these intervals 
by their median redshifts, z = 0.45, z = 0.55, and z - 0.75. 
Figure|2]shows the redshift distribution of our sample of galax- 
ies. 

Dust reddened star-forming galaxies are estir nated to make 
up ~ 2 0% of galaxie s in low-redshift c luster s ( Strateva et al 
200 ll: iGavazzi et all 120021; iBell et all 12004 iFranzetti et si 



0.45 



0.55 



0.75 




■ 1 



Spectroscopic Redshift 



Fig. 2. Redshift distribution of the red-sequence galaxies in the 
3 redshifts-bins with no or negligible emission lines, before any 
signal-to-noise selection. 



20071) . This contamination is likely to be more important in 
the strong-emission line regime To avoid contamination from 
these galaxies that do not have red stellar populations, and due 
to the difficulty of measuring reliable absorption line indices in 
galaxies with emission lines, we discard, for our absorption line 
analysis, galaxies showing evidence of emission lines. Our se- 
lection is based both on visual inspection of the 2D- spectra and 
on measurements of the line equi valent widths (EW s). These 
measurements were taken from Po ggianti et al. I (l2006h . We also 
measured the [Oil] and H/? emission-line EWs after carefully 
subtracting the underlying stellar continuum (see Moustakas & 
Kennicutt 2006 for details). After comparing the visual inspec- 
tion with the quantitative measurements, we keep galaxies with 
equivalent widths of [OII]/l3727 < 7 A. In Fig.[T]we distinguish 
between galaxies with strong emission lines and those with no- 
or weak- emission lines (N-i-W, hereafter) in the spectroscopic 
sample of red-sequence galaxies. Figure [3] shows the fraction 
of N-hW red-sequence galaxies as a function of redshift and 
structure velocity dispersion. The scatter is large and there is 
no trend with either crduster or redshift. On average, 24% of the 
red-sequence galaxies show EW[OII]3727>7A. Section[3]dis- 
cusses the nature of the emission in these galaxies. 

We also remove the brightest cluster galaxies (BCGs) from 
the sample, as their position in the cluster may lead to a differ- 
ent evolution from the rest of the red population (De Lucia & 
Blaizot 2007). A photometric study of the BCGs in the EDisCS 
clusters is presented in Whiley et al. (2008). 

Our initial sample includes a total of 337 N-hW red- 
sequence galaxies. Their mean spectral signal-to-noise ratio per 
A (S/N), measured between 4000 and 4500 A (rest-frame) 
is ~ 17, 19 and 12 per A at redshift 0.45, 0.55 and 0.75, re- 
spectively. We consider all galaxies independently of their dis- 
tance to the cluster centers except for the morphological anal- 
ysis. The radial cuts, when applied, will be stated explicitly. 
However, to measure line-strength indices, we require a mini- 
mum signal-to-noise ratio of 10 per A , measured between 4000 
and 4500A (rest-frame). Thie leave us, for our spectroscopic 
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Fig. 1. V, I colour-magnitude diagrams for the 24 EDisCS structures considered in this study. Galaxies with no- or weak- emission 
are represented by red cycles, while galaxies with strong emission are marked as blue circles. The fit, with a fixed slope of -0.09, 
of the secure cluster members with no or only weak emission is shown as a solid line. The dashed lines trace the ± 0.3 mag 
boundary taken into account to select the red-sequence galaxies. This corresponds to the measured dispersions of cll232. 5-1250 
and cll216. 8-1201, our two most populous galaxy clusters. When photo-z membership is available, we show the photo-z selected 
cluster galaxies as small black dots. 
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Redshift cr(km/s) 

Fig. 3. Fraction of red-sequence galaxies that have no or only weak emission lines (EW[OII] <7A) as a function of redshift and 
group/cluster velocity dispersion. This fraction is calculated as the ratio between the N+W red-sequence galaxies and the total 
number of red-sequence galaxies in the spectroscopic sample. 



analysis, with a total of 215 galaxies satisfying our selection 
criteria, i.e, 70% of the initial sample (see Table[T]for their dis- 
tribution). 



3. Source of ionization for the emission line 
red-sequence galaxies 

As inferred from Fig. [3l a non negligible fraction of red- 
sequence galaxies exhibit emission lines, with equivalent 
widths reaching ~30A for [Oil] and ~12A for H/3. Although 
these galaxies are discarded from further investigation, we are 
interested in trying to unveil their nature. We investigate below 
the various possible ionization sources and look for evidence 
of evolution with redshift. 

We perforrn a similar analysis than lYan et al. (l2006h . 
Yan et alJ (1200 6) studied a large sample of red galaxies from the 
Sloan Digital Sky Survey (SDSS) finding emission in 52.2% of 
them. They showed that they could be classified in two main 
groups: galaxies with high [OII]^3727A/Ha (LINER-like and 
quiescent galaxies, representing 20.6% of the total red populu- 
tion with emission) and low-[OII]/Ha galaxies (mostly dusty- 
starforming ones, with a small fraction of Seyferts, accounting 
for 9% of the red-sequence galaxies with emission). Counting 
galaxies with only [Oil] seen in emission (instead of both [Oil] 
and Ha) increases the fraction of LINERS to 28.8%. Finally, 
14% of the red galaxies have Ha detected but no [Oil], mak- 
ing them difficult to classify. Nevertheless, Yan et al. (2006) 
show that the majority of these galaxies have [NII]/Ha > 0.6 
suggesting a non-starforming origin for their emission lines. 

In the present case, the determination of the ionization 
sources is hampered by the fact that the wavelength cover- 
age of the spectra does not include Ha. We consider H/3 in- 
stead. For each galaxy, the stellar continuum has been fitted 
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Fig. 4. Ratio between [Oil] and H/3 for the full EDisCS sam- 
ple of cluster members (black dots). Red, green, and blue stars 
indicate the red-sequence galaxies in the 0.45, 0.55, and 0.75 
redshift bins, respectively. The[OII]/HyS distribution is also pre- 
sented in the form of a histogram: the black-filled zone rep- 
resents the red-sequence galaxies, to be compared to the full 
sample. Dashed lines indicate the [OII]/Hj6 = 6.7 limit above 
which the galaxy emission is considered to be due to mecha- 
nisms other than star formation. 



and the [Oil] and H/3 equivalent widths subsequently measured. 
FigurelHillustrates the distribution of [Oil] and H/3 EWs for the 
red-sequence galaxies and for the full EDisCS spectroscopic 
sample. 
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Out of the EDisCS full sample of red-sequence galaxies 
(393), 243 have spectra with a wavelength range appropriately 
centered and wide enough to include both [OII]/l/13727A and 
Hy6. The redshift distribution of these galaxies has a mean of 
0.52 and a standard deviation of 0.09. We count 57 gala xies 
with a 2-cr detection (threshold identical tol Yan et aL ( 20061) ) of 
both [Oil] and H/3, 68 with only one of the two lines detected, 
and 118 wi th no detection at all (at the 2-cr level). Hence, very 
similarly to lYan et al. I (l2006h . we obtain 5 1 .4% of emission line 
red-sequence galaxies. 



Following iKewlev et al.l (120041) . we consider that star for- 
mation cannot induce [OII]/Ha larger than 1.5. This value is 
derived before any reddening correction and wa s independently 
confirmed by iMoustakas & Kennicutti (120061) . We assume a 



mean value of Ha'/H/?=4.46, which was derived by lYan et al 
t2006.) before any reddening correction, and which corresponds 
to a median extinction of Ay =1.40, assuming Type B recombi- 
nation. Under this condition, [OII]/HjS= 6.7 draws a boundary 
between star forming galaxies and those w ith emissio n pow- 
ered by an AGN. As shown in figure 11 of lYanet all f2006). 
this limit is very conservative. Indeed, one hardly finds any star 
forming galaxies with [OII]/H;6 > 5. However, it does elim- 
inate most Seyferts and selects mostly LINER-type galaxies. 
Interestingly, the limit of 6.7 seems to represent a natural upper 
limit for the bulk of the whole cluster galaxy population (see 
the histogram part of Figure |4|l. For our statistics below, we 



use the empirical bimodal demarcation advised by lYan et al 
(|2006») . i.e., EW([OII])=18EW(H/3)-6 rather than a constant 
[OII]/I^ ratio. We have checked, however, that both methods 
lead to identical results (within 1-2%). 

Table |2] shows the emission-line properties of the red- 
sequence galaxies with spectra wide enough to cover li/3 and 
[Oil]. To go beyond upper and lower limits we would need, 
at least, another set of two lines to build more diagnostic di- 
agrams (e.g., Kewley et al. 2006). However, we believe that 
our counts give reliable hints on the major trends for the ion- 
ization processes at play in our sample. The fraction of qui- 
escent red galaxies is very m uch the same as the percentage 
(120061) in the local Universe. The obvi- 



reported by lYan et al. 
ous difference between the two studies resides in the nature 
of the galaxies which d o have detections of both H/3 and [Oil] . 
While Yan et al.l ( 2006b find that only 9% of those galaxies have 
low-[OII]/HjS ratios, we obtain more than twice this percent- 
age, 19%. In addition, only 4% of our red-sequence galaxies 
have high- [OII]/HjS ratios characteristic of LINERS. This frac- 
tion might increases up to 8% if we include galaxies with [Oil] 
detected and H/S undetected. However, this is still a factor of 3 
smaller than the fracti on of high-[OII ]/H6 galaxies in the red- 
sequence reported by Yan et al. ( 20061) . 29%. 

The conclusion for this tentative identification of the ion- 
izing sources in the EDisCS emission line red-sequence galax- 
ies is that most of them are dusty star forming ones. We can 
not disentangle, at this stage, whether time or environment, or 
even more technically the use of li/3 instead of Ha, leads to the 



[Oil] 




[0II]/I^8 


Fraction 


N 


N 




48.5 % 


N 


Y 


uncertain 


24% 


Y 


N 


high 


4% 


Y 


Y 


high 


4% 


Y 


Y 


low 


19.5% 



Table 2. Emission line properties of the red-sequence galax- 
ies for which the spectra covers both [Oil] and H/3. The 
first two columns indicate whether or not [Oil] and H/3 have 
been detected (at 2-cr level). High and low [OII]/Hy6 stand for 
EW([OII])>18EW(Hje)-6and EW([OII])<=18EW(H/3)-6,re- 
spectively. 



4. Stellar velocity dispersions 

Velocity dispersions were measured i n all our spectra using 
the ID L routine PPXF developed by /Cappellari & Emsellem 
(120041) . This routine, based on maximum penalized likelihood 
with optimal template, is specially suited to extract as much 
information as possible from the spectra while suppressing 
the noise in the solution and, therefore, is perfect to measure 
the kinematics in low S/N data. This algorithm estimates the 
best fit to a galaxy spectrum by combining stellar templates 
that are convolved with the appropiate mean galaxy velocity 
and velocity dispersion. The final values of these parameters 
are sensitive to the template missmatch and, therefore, the use 
of this technique requires templates which match closely the 
galaxy spectrum under scrutiny. This is achieved with the use 
of extensive stellar library spanning a large range of metallic- 
ities and ages. We use 35 synthetic spectra from the library 
of single stellar population models from Vazdekis et al. (2009, 
in prepa ration), which makes use of th e new stellar library 
MILES dSanchez-Blazquez et alj|2006a) . degraded to the res- 
olution of EDisCS spectra. The library contains spectra span- 
ning an age range from 0.13 to 17 Gyr and metallicities from 
[Z/H] = -0.68 to [Z/H] = -1-0.2. Errors were calculated by 
means of Monte Carlo simulations, in which each point was 
perturbed with the typical observed error, following a Gaussian 
distribution. Because the template mistmach affects the mea- 
sure of the velocity and cr determined with PPXF, a new opti- 
mal template was derived in each simulation. The errors were 
obtained as the standard deviation of a total of 50 simulations. 
The mean uncertainty in the velocity dispersion calculated in 
our spectra with S/N > 10 is ~25%. Given the signal-to-noise 
ratio and the resolution of our spectra, cr values below 100 
kms ' are considered untrustworthy and, consequently, galax- 
ies with cr < 100 kms ' are eliminated in all analysis involving 
the comparison with cr. 



5. Line-Strength indices 

For all of the galaxies in the re d sequence, we me asure 



difference with the results of I Yan et al.l (|2006). Further inves 



tigation, both in field intermediate redshift galaxies and in low 
redshift cluster galaxies, is required to shed light on this matter 



Lick/IDS indices as defin ed by iTrager et alJ (119981) and 
Worthey & Ottavianil (Il997h . The strength of the Lick/IDS as 
a function of age, metallicity, and individual chemical abun- 
dances has been calibrated by many authors for a simple stellar 
population, i.e. a population formed in a single, instantaneous 
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burst (e.g. lWorthevIll 994 IVazdekisIll 9991: Fhomas et alJlioOS 
TMB03 here after), and for more com plicated star formation 
histories (e.g. iBruzual & Charlotil200 3). These indices remain 
the most popular way to extract information about the stellar 
ages and metallicities from the integrated light of galaxies. 

The errors on the indices were estimated from the un- 
certainties c a used by photon noise using the formula by 
(119981) and from uncertainties due to the wave- 



Cardiel et al 



length calibration, derived using Monte-Carlo simulations. The 
variances were estimated from the residuals between the ob- 
served spectrum and the "best template", which was obtained 
in the calculation of the velocity dispersion, and which was 
previously broadened to match the line width of each galaxy. 
For the clusters with z > 0.6, indices redward of 4500 A were 
affected by telluric absorption in the atmosphere and by sky 
subtraction residuals and were not measured. For the clusters 
with z ~ 0.45 all Lick/IDS indices from HSa to Hfi were mea- 
sured, as well as D4000. Each individual spectrum was visually 
examined to check for indices affected by sky residuals or tel- 
luric absorption in the atmosphere. All affected indices were 
discarded from any subsequent analysis. 

Line-strength indices are sensitive to the line broadening 
due to both instrumental resolution and the stellar velocity dis- 
persion. In order to compare the galaxy spectra with stellar 
population models and to compare the line-strength indices of 
galaxies with different velocity dispersion, the indices need to 
be corrected to identical levels of intrinsic Doppler and instru- 
mental broadening. 

Because the models we are using (see Sec. |6]l predict, not 
only line-strength indices, but the whole spectral distribution, 
we can degrade the synthetic spectra to the resolution of the 
data. We decided to broaden all of the observed and synthetic 
spectra to a final resolution of 325 km s"' (including the ve- 
locity dispersion of the galaxy and the instrumental resolution) 
before measuring the indices. Galaxies with a velocity disper- 
sion higher than ~ 315 km s ' could not be broadened, but the 
number of galaxies with velocity dispersions above this limit 
is very small and, therefore, we decided to include them on the 
plots, although they have not been included in the data analysis. 

One of the problems of using Lick/IDS-based models (e.g., 
Worthey 1994; Thomas, Maraston & Bender 2003) is that the 
data need to be transformed to the spectrophotometric system 
of the Lick/IDS stellar library. To do this, it is common to ob- 
serve stars from the Lick/IDS library using the same instru- 
mental configuration than for the science objects and to derive 
small offsets between the indices measured in those stars and 
the ones from Lick/IDS. When analyzing data at high redshift, 
this is of course impossible. In principle, no further corrections 
to the indices are required when comparing data with stellar 
population models using flux-calibrated libraries. Although we 
do not have to apply the offsets to our indices, we have com- 
puted them using the stars in common between MILES and the 
Lick/IDS library because it may be useful for other studies. The 
final offsets and the comparison can be found in Appendix lApl 



Some of the Lick/IDS indices are affected by emission 
lines. In particular, emission, when present, fills the Balmer 
lines, lowering the values of the indices and, hence, increas- 
ing the derived age. High-order Balmer indices (H6 and Hy) 
are much less affected by e mission than the classic index H/3 
( Worthey & Ottavianill 199 71). but nevertheless still affected. To 
correct for emission, two different approaches are normally 
adopted: First, assume a correlation between the equivalent 
width of some other emission line and the emission in the 



Balmer line (e.g. iTrager et al.ll2000ah . Second, fit an optimal 
template and subtr act the emission directly from the residuals. 
Nelan et alJ (120051) have shown that the first approach suffers 
from significant uncertainties, in part because there are almost 
certainly several competing sources of ionization in early type 
galaxies. The second approach requires that the fit of the opti- 
mal template to the individual spectrum is very good. The S/N 
of our individual spectra simply do not allow us to explore that 
option. Therefore, instead of trying an emission correction to 
the Balmer lines, we have analyzed the differences in the results 
by first eliminating all the galaxies showing any emission. It is 
re-assuring to see that none of our conclusions change when 
the (weak) emission line galaxies are excluded (we remind the 
reader than the galaxies with strong emission lines have been 
excluded from the analysis). 

Table [3] available at the CDS, lists the measured line- 
strength indices in our sample of galaxies with non-or weak 
emission lines. A portion of the table is shown here to show its 
content an structure. 

5.1. The local sample & aperture correction 

To increase the time baseline of our analysis, we com- 
pare the EDisCS spectra with 36 early-type galaxies (ellip- 
ticals and lenticulars) in the Coma cluster at redshift z - 
0.02. The characteristi c s of th is sample are described in 
Sanchez-Blazquez et al. (l2006bl) . Because the galaxies were 
selected morphologically we checked first that all the objects 
belong to the red-sequence of this clu ster, u sing the colours 



and magnitudes from Mobas her et al. We are aware 



that the comparison with the Coma cluster is probably not the 
most appropriate, as the the velocity dispersion of this clus- 
ter exceeds that of all the clusters in our sample at interme- 
diate redshifts. For this reason, we will base all the conclu- 
sions of this paper in the intercomparison of the EDisCS clus- 
ters. The local sample is included here to show that some of 
our results can be extrapolated to the local Universe, where 
higher S/N spectra can be obtained. One of the main prob- 
lems when comparing observations at different redshifts is that, 
for a fixed aperture, one is sampling different physical regions 
within the galaxies. Early-type galaxies show variations in their 
main spectral characteristics with radius and, therefore, a cor- 
rection due to these aperture differences is necessary. In or- 
der to compare directly with the sample at medium- and high- 
redshift, we extracted the ID spectra in the same way as the 
EDisCS spectra, inside an aperture equal to the FWHM of the 



^ Note that these offsets will not correct for any systematic effect 
due to a not-perfect calibration in the data. They are only useful to cor- 



rect the fitting-function-based models from the non-perfect calibration 
of the Lick/IDS stars. 
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spatial profile. However, aperture effects are not entirely miti- 
gated via such an extraction, as the EDisCS spectra were ob- 
served with a l"-wide slit, equivalent to a much larger physi- 



cal ap er ture at the redshif t of t he Coma cluster J 0rgen s en et al 



I995I) . |j0rgensenl (Il997h and |j0rgensen etal] (l2005l) derived 



aperture corrections using mean gradients obtained in the lit- 
erature. The index corrected for aperture effects can be ob- 
tained as log(index)coiT = log(index)ap -1- a log ^ for those 
indices measured in A and indeXcorr = indeXap + a log ^ 
for molecular indices measured in magnitudes, where r„or is 
the final normalised aperture, and rap is the equivalent cir- 
cular aperture, obtained as 2 rap - 1.025 x 2(xy/;7r)'^^, be- 
ing X and y the width and length of the rectangular aperture. 
J0rgensen (1995, 1997) calculated a parameters for a large sub- 
set of Lick indices but, in some cases, using the gradients mea- 
sured in a very small sample of galaxies. In this work, we have 
taken advantage of the large sample o f galaxies with measured 



line-str en gth gradients published in [ Sanchez- Blazque z et al 
( 2006al) . ISanchez-Blazquez et al.l ( 120071) and Jablonka et al 



J2007h to calculate new aperture corrections. These also in- 



clude corrections for the higher-order Balmer lines that have 
not been published before, as far as we are aware. AppendixlB] 
lists the new a parameter calculated in this paper. We correct 
all our indices to mimic the physical aperture of the slit used 
to observe the galaxies at z = 0.75. The aperture corrections 
(a log — ) for all the indices at the different redshift bins are 
Usted in"Table|4l 

Line-strength gradients show la rge variations be- 
tween early-type galaxies (see, e.g. Carollo et al. ' '19931 
iGorga s et al. 1990; D a vies et al. 1993; Halliday 1998 ; 
Sanchez-Blazquez et al. I l2006at iKuntschner et alj l20o5 
Jablonka et al.l 120071 among others), and the slope of the 
gradients does not seem to coiTelate clearly with any of their 
other properties. For this reason, we use the mean gradient to 
compute the aperture coiTections for all the galaxies. However, 
in order to calculate the error in the aperture coiTection due to 
the fact that early-type galaxies show a large variation in the 
slope of their gradients, we perform a series of Monte Carlo 
simulations where the mean index-gradient was perturbed by a 
random amount, given by a Gaussian distribution with cr equal 
to the typical RMS- dispersion in the gradients The errors in 
the aperture correction calculated this way are also indicated in 
TablelH In some cases (e.g., higher-order Balmer lines), these 
errors are very large. In fact, despite the large mean aperture 
correction for H6 and Hy indices, the aperture correction is 
compatible, within those errors, with being null. These errors 
have been added quadratically to the eiTors of the indices for 
the Coma galaxies. For the 0.45 and 0.55 bins, the error in 
the aperture correction due to the scatter among the mean 
gradients is negligible (less than 1% of the total coiTection) 
and we do not list them. 

Gradients might evolve with redshift. That is the reason 
why we coiTect the indices in the Coma cluster to the same 
equivalent aperture as the highest redshift bin instead of cor- 
recting the highest redshift measurements as is commonly 
done. This way, the correction is performed to the galaxies 



ml m 



Fig. 5. Velocity dispersion distribution of the individual spec- 
tra stacked in the different redshift bins. From bottom to top: 
z=0,0.45,0.55,0.75. 



where the gradients used to derive the aperture correction have 
been measured. 



5.2. Stacked spectra 

The degeneracy between age and metallicity effects on galaxy 
colours could well affect the interpretation of their colour- 
magnitude diagrams. This degeneracy can be partially broken 
by using a combination of two or more spectral indices chosen 
for their different sensitivities to age and metallicity. The S/N 
of our galaxy spectra is not high enough to do this analysis on 
an individual basis. Therefore, in order to derive the evolution 
with time of the galaxy mean ages and chemical abundances, 
we stacked the galaxy spectra in each of the 0.45, 0.55 and 0.75 
redshift bins. 

We made sure that the distribution of galaxy velocity dis- 
persions (cr) in all the different redshift intervals is similar 
(Fig. |5]). This a prerequisite to any further comparison be- 
tween the three redshift bins, as most line-str ength indices ar e 
strongly correlated with this parameter (e.g. iGonzalej 1993 ). 
To explore the possible dependence of star formation history 
on galaxy mass, we distinguish between galaxies with cr higher 
and lower than 175km s ' in all three redshift groups. We per- 
form a Kolmorogov-Smirnov test comparing the distribution 
of at different redshifts to look for possible significant dif- 
ferences. We performed the test separately for galaxies with cr 
larger and smaller than 175 km s ' and found them compatible 
with being drawn from the same distribution. 

Before adding them, we normalize all spectra by their mean 
flux in the region between 3900A and 4400A, to avoid a bias 
towards the indices of the most luminous galaxies. Finally, we 
co-add them, clipping out all the pixels deviating more than 2-cr 
from the mean. Figure |6] displays the resulting 6 stacked spec- 
tra for the three redshift bins and the two velocity dispersion 
regimes. 
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Fig. 6. Final stacked spectra for our three redshift and velocity 
dispersion bins. The shaded areas identify the central bands of 
the four indices used in our spectral analysis. 



Earlier works have studied the evolution of red-sequence 
galaxies as a function of magnitude (e.g. De Lucia et al. 2004; 
De Lucia et al. 2007; Rudnick et al. 2008, in preparation). 
These studies have put the boundary between faint and bright 
galaxies at Mv= -20. Assuming a formation redshift of 3 fol- 
lowed by passive evolution, our cr=175 km s"' cut corresponds 

to Mv 20.0, -20.2 and -20.2 mag at z = 0.45, 0.55 and 

0.75, respectively (see Figure |7]). As can be seen, the magni- 
tude cut is not very different between our study and the above 
mentioned photometric studies. However, those studies reach 
much fainter magnitudes in their "faint" bin. 

Velocity dispersions and line-strength indices for the 
stacked spectra were measured using the same techniques as 
for the individual spect ra. Errors on the in dices were calcu- 
lated with the formula of lCardiel et alj (Il998h . using the signal- 
to-noise of the stacked spectra. These errors do not reflect 
the differences between the spectra added in each bin, which 
are larger than the formal errors obtained using the signal-to- 
noise ratio. However we do not intend to study the distribution 
of galaxy properties but rather their mean values. To explore 
the robustness of our results to errors in velocity dispersion, 
we performed 20 Monte Carlo simulations in which the in- 
dividual galaxy velocities were randomly perturbed following 
a Gaussian error distribution. After each simulation, the spec- 
tra were stacked again, indices were measured, and their mean 
and RMS -dispersion were calculated. Fig. [8] compares the in- 
dices obtained in the original stacked spectra with the mean in- 
dices obtained from the 20 simulations. The separation in low- 
and high-o" galaxies is robust to the errors in velocity disper- 
sion. When the RMS -dispersion for all the simulated galaxies 
is larger than the error we have calculated for the individual in- 
dices, we add the quadratic diff'erence as a residual error to the 
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Fig. 7. Relation between the absolute magnitude Mv and the 
velocity dispersion for our sample of galaxies in the three con- 
sidered redshift bins. The right label shows the absol ute mag- 
nitude Mv evolved passively to z = using Bruzual & Chariot 
(l2003h models at solar metallicity. The vertical dashed Une in- 
dicates the position corresponding to our cr-cut (175 km s"') 
while the horizontal line marks the mean magnitude for the 
galaxies with a sigma between 160 and 180 km 

original error in the index. We have also checked that the mean 
values of the individual indices for all the stacked galaxies is 
compatible, within the errors, with the index measured in the 
stacked spectra. 

The Lick/IDS indices measured in the stacked spectra as 
well as their errors are listed in Table |5] 



6. Stellar population models 

In this study we compare our measured values with the 
models by Vazdekis et al. (2008, in preparation, V08 here- 
afteJl). These models are based on the MILES stellar library 
(ISanchez-Blazquez et al. ll2006cV The spectra from this libiary 
have a resolution of 2.3A , constant along the whole wave- 
length range, and are flux-calibrated. Thanks to this, the models 
predict not only individual indices, but the whole spectral en- 
ergy distribution from 3500 to 7500A as a function of age and 
metallicity, for ages between 0.1 and 17.78 Gyr and metallici- 
ties from [Z/H]=-1.68 to -hO.2 dejfl 

However, the chemical composition of stars in the MILES 
library (as in the rest of empirical libraries included in stellar 
population models) matches that of the solar-neighborhoodH 
This means that the models, in principle, can only predict ac- 
curate ages and metallicities for stellar populations with the 



The 



models 



are 



publicly 



available 



at 



=0.02 



http://www.ucm.es/info/Astrof/miles/models/models.html 

* In the normalization to the solar values a solar metallicity Zq 
has been used 

^ Although the detailed chemical composition of all the stars in the 
MILES is unknown, it is reasonable to assume that this is the case. 
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Fig. 8. Comparison between the indices measured in the 
stacked spectra and the mean index calculated in the 20 stacked 
spectra obtained with Monte Carlo simulations as described in 
the text to analyze the robustness of our separation of galax- 
ies with o" due to errors in this parameter The errors in the 
y-direction represent the rms of all the 20 measured indices, 
while the errors in the jc-direction represent the errors in the 
index due to the S/N of the stacked spectra. Big and small sym- 
bols show the indices measured in the stack with cr larger than 
175 kms"' and smaller than cr < 175 kms"' respectively. The 
meaning of the different symbols is indicated in the inset. 



same chemical patterns as the solar vicinity. We know that this 
condition doe s not apply to bright nearby early-type galaxies 
(see review bv lWorthevI 1998^ 1. where Mg/Fe and maybe Ti/Fe, 
Na/Fe, N/Fe and C/Fe are enhanced with respect to the solar 
values. This may result i n systematics errors in the me asured 
ages and metallicities (see Worthevl[l998 : Thomas et al.l l2003). 
In order to overcome this problem, we calibrate the models for 
different chemical compositio ns u sing the method first intro- 



duced by lTantaloetalJ(ll998l) and iTrager et al.l (l2000ah . This 



method uses stellar atmospheres to characterize the variation 
of the different indices to relative changes of different chemi- 
cal species. We first interpolate the existing models to make a 



grid with 168 ages from 1 Gyr to 17.8 Gyr in steps of 0.1 Gyr 
and 209 metallicities from -1.68 to -hO.4 in steps of 0.1. Then, 
we modify each index according to their fractional response to 
variations of different elements using (from Trager et al.); 



A/ 
T 



]~[(1+/?0.3(X/)) 



K/H]/0.3 



- 1, 



(1) 



where RQ o,{Xj) is the response function for element / at 
{XilH} - H-0.3 dex. As in Trager et al. , we also consider that 
the fractional light contribution from each stellar component is 
53% from cool giant stars, 44% from turnoff stars and 4% from 
cool dwarfs stars. We use the calibrations by Korn et al.l (l2005h 
instead of the ones by lTripicco & Belli (119951 (used by Trager 
et al.) because they include the higher order Balmer lines Wd 
and Hy, and they are computed at other metallicities different 
from solar. Ho wever, the differen ces between the two studies 
are very small (iKorn et al.ll2005l) . Following the approach of 
several previous authors, we do not change all the elements 
individually but we assume than some of them are linked by 
nucleosynthesis and, therefore, change them in lock-step. We 
built models where N, Ne, Na, Mg, Si and S are enhanced 
by 0.3 dex with respect to Fe while C is enhanced by +0.15 
dex. The reason that C is only enhanced by +0.15 is that en- 
hancing C by +0.3 dex brings the C/O ratio very close to the 



values at which a carbon star is formed (see iHoudashelt et al 



20021: iKorn et al.ll2005L for a discussion). Ca and the Fe peak 
elements are depressed to keep the overall metallicity con- 
stant. @To obtain response functi ons for enhancement s different 



from +0.3 dex we, again, follow lTrager et al.l (1200 0b): We first 
calculate A[Fe/H] and A[E/Fe] using A[Fe/H]=-AA[E/Fe] = 
-Y^A[E/H], where E represents the enhanced elements. A is 
the response of the enhanced elements to changes in [Fe/H] at 
fixed [Z/H]. For the chosen enhanced model, A - 0.929 (Trager 
et al. 2000). Then, we scale exponentially the response func- 
tions by the appropriate element abundance. We built models 
with enhancements from [E/Fe]=-0.1 to 0.5 in steps of 0.05. 

Because we are not working a t the Lick/IDS reso lution, in 
principle, we should not be using iKorn et al. (2005) response 
functions to compute the sensitivity of different indices to vari- 
ations of chemical abunda nce ratios, as they were computed 



at this resolution. However. iKorn et al. 1(12005') showed that the 



influence of the resolution on the response functions is very 
small. 

Note that the models used here (as well as most exist- 
ing models in the literature) do not include isochrones with 
chemical abundance ratios different from solar. The effect of 
the different chemical abundance ratios is only included in the 
atmospheres. The influence that the inclusion of a-enhanced 
isochrones may have in the final predictions is still unclear, 
especially at super-solar metallicities, but some studies have 
shown that it is much smaller than the one in the model atmo- 



^ Despite the fact that Ca is an a-element, theoretically linked with 
the overabundant elements as Mg, measurem ents show that it may be 
depressed in elliptical galaxies (although see IProchaska et al] jlOOSl) 
for a different point of view). We follow Trager et al. 2000 and include 
this element in the depressed group. 
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spheres (ICoelho et al .l2007h . at least for line-strength indices at 
wavelengths shorter than ~5750. 

We want to caution the reader that the absolute values of 
ages obtained directly with stellar population models are sub- 
ject, not only to errors in the data (which are taken into account) 
but also to errors in the stellar population models. Stellar evo- 
lutionary models are still affected by un certain ties that leave 
room for improvement (see, e.g. ICassisil 2005 ). as illustrated 
by the non-negligible differences still existing among the re- 
sults provided by different theoretical groups. Furthermore, it 
must be noted that the oldest ages of the models we are us- 
ing (17.78 Gyr) are older than the current age of the Universe. 
The derivation of very old spectroscopic ages, inconsistent with 
the ages derived from colour-magnitude diagrams for globu- 
lar clu sters and with the ag e of the Universe, was first pointed 
out bv lGibson et al.l (Il999h and is a well-known problem in the 
community. Some solutions have been proposed, such as the in- 
clusion of atomic diffusion in the stell ar evolutionary mode ls or 
the use of a-enhanced isochrones (see Vazdekis et al.l2001 ), al- 
though none of these have bee n implemented in publ icly avail- 



able stellar population models. ISchiavon et al.l (12002 ) showed. 



for the particular case of 47Tuc, that the luminosity function of 
the red giant branch is underestimated in the stellar evolution- 
ary models and that the use of the observed luminosity func- 
tions instead of theoretical ones results in derived ages for this 
cluster consistent with the age of the Universe and with those 
derived directly from the colour-magnitude diagram. A similar 
effect in super-solar metallicity models would cause spectro- 
scopic ages to be 30% too high (see Schiavon et al. 2002 for 
details). It is beyond the scope of this paper to discuss these 
uncertainties in stellar population models, but it is clear that 
conclusions based on the absolute values of age are simply 
too preliminary. Relative differences in age and metallicity are 
more reliable than absolute values, and we will therefore base 
our interpretations on relative values. 

Lick/IDS indices, as well as colours, respond to variations 
in both age and chemical abundances, but the relative sensi- 
tivity of the different indices to these two parameters is not 
the same. In order to break the existing degeneracy between 
age and metallicity, it is necessary to combine two or more in- 
dices. In what follows, for the main analysis, we use 4 indices: 
Fe4383, CN2, HSa and Hyf. We choose these because, of all 
the Lick/IDS indices that could be measured in all galaxies, 
they are the most sensitive to variations of age (¥16 and Hy) and 
total metallicity (CN2, Fe4383). TMB03 recommend CN and 
Fe4383 as the best blue indices to calculate a/Fe abundance 
ratios. Furthermore, we choose Hjf instead of Hja because 
the ages measured with Hja are systematically younger than 
those measured with Hyp or ¥1/3 when a-enhanced models are 
used ( Thomas et alJl2004 ). As th e origin of these diffe rences is 
unclear, we follow the advice by Thomas et al. ( 2004 1) and use 
Hjf, despite the higher photon noise in this index compared to 
its wider version (Hy^)- 



lected indices (HSa, Hy/r, CN2, and Fe4383) with the models 
by V08 using a minimization routine. We find that the de- 
rived ages are independent of the chosen Balmer index (H(5a 
or Hyp)- The resulting ages, metallicities, and [E/Fe] are listed 
in Table |6] Index-index diagrams comparing the indices mea- 
sured in the stacked spectra with the predictions for single stel- 
lar population models by V08 are shown Fig.|9] In each panel, 
[E/Fe] is chosen to be as close as possible to the values given 
in Tabled 

Although the selection of our sample has been done us- 
ing colours, many previous studies of the red-sequence were 
restricted to morphologically classified early-type galaxies. To 
compare with those, the analysis in this section has been re- 
peated for the subset of red-sequence galaxies that are morpho- 
logically classified as E or SO galaxies. For the clusters imaged 
with the ACS (see Table [T]|, we use their visual classification 
dOesai et al. 2007b . For the others, we apply the selection cri- 
terion of Simard et al. (2008), based on the bulge fraction de- 
rived from our VLT 1-images (see Sec|9]for the criteria to select 
early-type galaxies). 

The last three rows of Table |6] show the SSP-parameters 
derived for the red-early-type galaxies. The first thing that can 
be seen is that the SSP-parameters for the whole red-sequence 
and for the morphologically classified early-type galaxies are 
compatible within the errors. Therefore, we will not discuss 
further the possible differences and will restrict our analysis to 
the complete red-sequence sample. 



7.1. Massive galaxies 

It can be seen in Table |6] and Fig. |9] that the age difference 
between the redshift bins of the most massive galaxies (cr > 
175 kms"') corresponds to the expected age difference of the 
Universe at those redshifts. The absolute age obtained from the 
models for these galaxies corresponds to a redshift of 1 .4. In 
other words, they are compatible with being formed at z > 1 .4 
and evolving passively since then. The metallicity measured 
from CN2 does not evolve either, as expected in a passive evo- 
lution scheme. 

When we look at the panels using Fe4383 instead of CN2 
we can see that the most massive galaxies at z = 0.75 seem 
to have a lower Fe43 83 than predicted by passive evolution. 
^J0rg ensen et al. (2005) similarly found a weaker Fe4383 index 
for more than half of the galaxies in a sample of red-sequence 
galaxies (including those with cr larger than 175 km s"') in the 
cluster RX J101522.7-1357 at redshift z = 0.83. This result, if 
confirmed, could be indicating that at least some massive galax- 
ies have experienced some chemical enrichment since z = 0.75. 
It will be necessary to measure other Fe-sensitive indices to 
confirm this trend. Unfortunately, the wavelength range cov- 
ered by our spectra does not allow us to do this and, therefore, 
we will not discuss this further. 



7. Age & Metallicity 

We derive age, metallicity and the ratio of a-element enhance- 
ment [E/Fe] for our stacked spectra comparing our four se- 



7.2. Less-massive galaxies 

Contrary to their massive counterparts, galaxies with cr < 
175 km s ' do not show any evolution, within the errors, in 
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Fe4383 Fe4383 

Fig. 9. Index-index diagrams for the indices measured in the stacked spectra at different redshifts. Left panels (big symbols) show 
measurements on the stacked spectra of galaxies with cr larger than 175 km s"' while right panels (small symbols) show the same 
measurement for galaxies with cr smaller than 175 km s The different redshifts bins are represented as different colours (and 
different symbol shapes); z - 0.75: blue triangles; z - 0.65: green stars; z - 0.45: red circles; z — 0.02: orange squares. Model 
grids of constant age (solid lines) and constant metallicity (dashed lines) from V08 with [E/Fe]=-i-0.35 (for the most massive 
galaxies ) and [E/Fe]= 4-0.1 (for the lower sigma bins) are superimposed. The ages and metallicities of the different models are 
indicated in the Figure. 



age or in metallicity between z - 0.75 and z - 0.45. They 
also evolve less than expected from a pure passive scenario 
between z - 0.45 and z. — .0. A similar result was re- 
ported by ISchiavon et al. I (l2006h with a field galaxy sample. 
Comparing galaxies from the DEEP2 survey around z - 0.8 
with an SDSS local galaxy sample, they showed that the W5f 
variation was less than predicted by passively evolving models. 
However, they do not group galaxies by mass. Along the lines 
discussed by the authors, our results suggest that either indi- 
vidual low-mass galaxies experience continuous low levels of 
star formation, or that the red-sequence is progressively built- 
up with new and younger small galaxies. This latte r hypothe 



sis is supported by a number of recent works (e.g. iBell et al 



20041 ; IPoggianti et al.ll2006HDe Lucia et"ai]|2007l; iHarker et al 



20061; iFaber et aLll2007h . This differential evolution for mas- 
sive and less massive galaxies implies that the mean differ- 
ence in luminosity-weighted mean-age between massive and 
intermediate-mass red-sequence galaxies increases with time. 



This should be taken into account when studying the evolution 
of the colour-magnitude relation with redshift. 

Low-cr galaxies show a lower value of [E/Fe] than high-cr 
galaxies, in agreement with the results at low redshift. Most of 
the "enhanced" (E) elements we are considering here arise from 
massive stars and the bulk of their mass is released in Type II 
supernova explosions, while at least 2/3 of the Fe is released to 
the interstellar medium by Type la supernovae, with a delay of 
~ 1 Gyr Therefore, the ratio [E/Fe] has classically been used 
as a cosmic clock to measure the duration of the star formation. 
The lower [E/Fe] in low-cr galaxies can be interpreted as a more 
extended star formation history, just as found locally. 

Table|6]also shows that, within the errors, [Z/H] and [E/Fe] 
do not vary with redshift for low-cr galaxies. Whether galax- 
ies have experienced low-level star formation or new galaxies 
have entered the red-sequence from the blue cloud, it is still un- 
clear how galaxies that have been forming stars until recently 
have the same chemical composition as those where star for- 
mation was quenched 2 Gyr previously (the lookback time be- 
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tween z~0.45 and z~0.75). Answering this question is not triv- 
ial, as there are multiple paths by which a galaxy may enter 
and leave the red-sequence. Furthermore, one has to be care- 
ful when interpreting single-stellar population equivalent pa- 
rameters. While galaxy mean ages are biased towards their last 
episode of star formation, the chemical composition depends 
more strongly on the oldest population (see Serra & Trager 
2007 for a quantitative study). Therefore, galaxies with differ- 
ent star formation histories can have similar [E/Fe] and [Z/H] if 
they have formed the bulk of their stars on similar time-scales 
and with a similar efficiency. 

8. Evolution of the index-cr relations with redshift 

The outcome of the previous section seems to contradict ear- 
lier studies, in which the evolution of the index-cr relation was 
found compatible with a high-redshift formation and subse- 
quent passive evolution for all red-sequence galaxies, irrespec- 
tive of their mass range (Kelson et al. 2001). We will now 
demonstrate that our index-cr relation is also consistent with 
this scenario, but that this analysis cannot provide robust con- 
straints on the star-formation histories of red galaxies. In fact, 
the evolution of the index-cr relation is also compatible with 
a more extended star formation history for the intermediate- 
mass galaxies if the progenitor-bias (i.e. the fact that the galax- 
ies with more extended star formation history (or quenched 
at later times) drop out of the red-sequence at high redshift 
van Dokkum et al. is taken into account. 

Figure [TO] displays the relation between Hja, Hjf, CN2, 
and Fe4383 and the galaxy velocity dispersion, for the Coma 
cluster and ou r EDisCS samp le 



Following iKelsonetaL] (1200 lb and lJ0rgensen et alj (120051) . 
among others, we first calculate the best linear fit to the data 
using the galaxy sample at z = 0.45. Then, keeping the slope of 
the relation fixed, we evolve its zero point with look-back-time. 
This approach assumes that, if galaxies were coeval and evolv- 
ing passively, the slope of the relation would not change. In 
reality, ;/ all the galaxies in the red-sequence were coeval and 
evolving passively, we would expect a small variation in the 
slope of the index-cr relation, as the variation of the indices is 
not completely linear with age. However, this variation would 
be very small over the range of ages considered. To calculate 
the lines of passive evolution, we used the V08 models, assign- 
ing solar metallicity to the galaxies with cr = 300 km s"', to 
match the observed value of the metal-sensitive index Fe4383 
at z = 0.45. We show the best linear fits in Fig.[TOl obtained by 
minimizing the residuals in the y-direction. We also show the 
expected evolution of the index-cr relation assuming three dif- 
ferent formation redshifts: Zf = 1.4, 2, and 3. The value Zf - 1.4 
is chosen as it corresponds to the mean age measured in the 
stacked spectra of the massive galaxies in the 0.75 galaxy group 
(see Sec. 15.2b . We only plot models with solar-scaled chemical 
abundances. 

Figure [TO] shows that a scenario where all stars form at 
Zf > 1 .4 and evolve passively afterwards is compatible with the 
observations. The exception is CN2 which we discuss below. 
In order to analyze this more quantitatively, we performed a 
f-test, comparing the linear fit values at cr = 200 km s ' (which 



is approximately the mean of the distribution of cr) with the 
ones predicted by passive evolution. The results are shown in 
Table|2l f-values higher than 1.9 would indicate that the proba- 
bility that we have rejected the null-hypothesis (prediction and 
measurement are equal) by chance is less than ~5%, i.e., a t- 
value higher than 1.9 implies that the passive scenario does not 
reproduce our relations within the errors. We do not obtain such 
f- values, implying that indeed a passive evolution model can re- 
produce the evolution of the zero-point of the index-cr relations 
from z = to, at least, z~0.75. Although we obtain a slightly 
better agreement between model and data if we assume a for- 
mation redshift above 2, we do not find any statistically signif- 
icant difference from a formation redshift of 1 .4. Noticeably, 
formation redshifts lower that 1 .4 produce an evolution of the 
index-cr relation significantly larger than the one observed. 

The relation of the other indices with cr is presented in 
Appendix |C] The analysis of these indices gives consistent re- 
sults without adding any new information to our study, but 
we show it for comparison purposes with eventual follow- 
ing studies. We did not include CN2 in the f-test procedure 
as the stellar population models with solar-scaled chemical 
abundances cannot reproduce the absolute values of this in- 
dex j_ThisJs__a_well known effect in nearby early type galax 



les ('Worthev'l99i ISanchez-BlazQuez et al.*2003l:lKelson et al 



2006; Schiavon et al 



2006; Graves et al. 2007) and it is at- 



tributed to differences in the chemical composition of these 
galaxies compared to the solar values. (Note that in Sect. 15.21 
we used models where C/Fe and N/Fe are enhanced with re- 
spect to the solar values in order to reproduce this index). 

Summarizing, we show here that, in agreement with previ- 
ous works, the evolution of the index-cr relation is compatible 
with a scenario where the red galaxies formed their stars at 
z>1.4 and have evolved passively since. However, this conclu- 
sion has been reached assuming that all red galaxies formed 
at the same time independently of their mass. The results 
of Sec. | 5. 21 argu e against this . Furthermor e, numerous works 
dOe Lucia et al. ! 2004, 2007; Kodama et a"l] |2004) indicate that 
the red sequence was not yet fully in place at z~l and that it 
has been growing since then. If this effect is taken into account, 
a more complex star formation history is allowed while keep- 
ing the slope and the tightness of the index-cr relations. This 
was proposed by Ivan Dokkum & Franxl (1200 lb to explain the 
evolution of the magnitude and the constancy and scatter of the 
col ours on the red sequ ence. 



Harker et al.l (l2006h showed that quenched models, where 



a constant star formation is truncated at evenly spaced time in- 
tervals can explain the evolution of the mean H6 of all the field 
red galaxies between z - 1 to z = 0. However, they did not 
explore the dependence of this evolution on galaxy mass. In or- 
der to study this aspect, we have used the iBruzual & Chariot 
(2003) population synthesis models with a Salpeter IMF and 
solar metallicity. We built a series of star-formation histories 
starting at Zf=2 where a constant rate of star formation of 1 
Mo/yr is quenched at different look-back times. Whenever a 
galaxy satisfies our criteria to belong to the red-sequence, we 
measure its indices. 

Figure [TT] shows the resulting line-strength indices assum- 
ing 3 different metalUcities, Z =0.001, 0.004 and 0.02. We 
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Fig. 10. Relation between line-strength indices and the velocity dispersion for the galaxies in the three different bins. Solid 
coloured lines represent the result of a linear fit to the data forcing the slope to be the same as the one obtained for the galaxies in 
the redshift bin at z = 0.45. The black lines represent the expected relation assuming that galaxies formed at Zf = 1 .4 (dotted-dashed 
line), Zf = 2 (thick, dashed line), and Zf =3 (dotted line), and have evolved passively since then and assigning a solar metallicity 
for galaxies with velocity dispersion of 300 km s"'. As can be seen, the evolution of the index-cr relation is compatible with 
a scenario in which the stars formed at z > 1 .4 and evolved passively since then, except for the CN2 index, which cannot be 
reproduced by any of the proposed models (see text for details) 

. However, the evolution of these relations are also compatible with more complex scenarios (see Sec.llOll. 



display, again, the relation of the indices with cr for each red- 
shift bin discussed in Sec. [8] While the index-cr relations of 
the most massive galaxies are only well reproduced by a sce- 
nario where star formation is truncated after 1 Gyr (i.e. the rest 
of the star formation histories produce relations systematically 
shifted compared to the observed one), less massive ones can be 
described by a variety of star formation histories, in agreement 
with our results of Sec. 15.21 Furthermore, this variety of star 
formation histories is supported by the larger scatter in metal- 
licity and age values at a given mass of low -cr galaxies com- 



pared to that of the most massive ones (e.g. Conc annon et al 



2OOOI) . Some studies have claimed that star formation histories 
of this nature would have problems reproducing the colours of 
red-sequence galaxies. As a sanity check (see Fig. [12]) we ver- 



ify that these star formation histories reproduce the colours of 
the red-sequence. 

This, of course, does not prove that different galaxies have 
different star formation histories, but it does demonstrate that 
more complicated scenarios than pure passive evolution are 
compatible with the index-cr relations. 



9. Morphological content of the red-sequence 

Table [T] lists the EDisCS structures wich have been observed 
with the HST/ACS in the F814W band, i.e, all six in the 0.75 
redshift bin, six at z = 0.55, and three at z = 0.45. In the fol - 
lowing, we use the visual classifications of iDesai et alj (l2007h . 
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Fig. 11. Evolution of the relation of line-strength indices with cr. We show with coloured, solid lines, the best linear fits to the 
observations, fixing the slope to that obtained for the redshift bin z = 0.45 (see Sec. [8]). Squares, diamonds, and black stars 
represent the predicted indices for a constant star formation rate starting at Zf = 2, for stellar populations with metallicities 
Z=0. 001, 0.004 and 0.02, respectively. At each of these metallicities, the sta r formation is truncated after 1 (largest symbols), 2, 
3, 4 and 5 (smallest symbols) Gyr. 



which has been performed down to /fl,„JZl=23 mag. Since we 
concentrate on the EDisCS spectroscopic sample, the actual 
magnitude cuts for the morphological classification are 1 (r < 
1") ~ 22 mag at z = 0.45 and l(r < 1") ~ 23 mag at z=0.75. 
This translates into restframe magnitudes of Mv~ -18.7 and 
~ -19.3 mag at z=0.45 and z=0.75, respectively. Once we ac- 
count for the dimming due to passive evolution (0.5 mag from 
z=0.75 to 0.45 assuming an age of 3.1 Gyr at z=0.75), the in- 
trinsic magnitude cut is the same at all redshifts. We note that 
some secondary structures at low redshift were discovered in 
high-redshift targeted fields. Consequently, they were observed 
down to lower magnitude limits than the primary targets at sim- 
ilar redshifts. As such, they are exceptions to the above general 
rule. However, they represent a small percentage of the total 
number of galaxies. 



Interestingly, only 15 E/SO galaxies out of 168 (i.e., 9%) 
are not on the red-sequence. This is similar to t he fraction of 



SExtractor Kron masnitudes 



blue early-type galaxies found at low redshift ( Bamford et al 
20081) . which is highly dependent on local density but varies 
between 12 and 2%. 

We concentrate on the N-hW red-sequence galaxies, the 
same for which we have analyzed the stellar population. 
Figure [T3] shows their Hubble type distribution. We also in- 
dicate the distribution of galaxies with cr < 175km/s. There is 
no apparent morphological segregation as a function of galaxy 
mass, i.e, whatever the galaxy mass range, the dominant mor- 
phological types on the red-sequence are clearly E and SO, al- 
though the entire Hubble sequence is covered. 

In order to look for possible evolution with time, we now 
restrict the sample to the EDisCS clusters (galaxy structures 
with cr larger than 400 km s"') and fix the area of analysis 
to the region falling inside a radius of O.6R200 from the clus- 
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Fig. 12. Measured V - / vs cr for the galaxies in the different redshift bins. Symbols are the same as in Fig.fTTl 
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Fig. 13. The distribution of Hubble types among the N+W 
emission line red-sequence galaxies in redshift bins: E(T=-5), 
S0(T=-2), Sa (T=l), Sab(T=2), Sd(T=7). The solid hne his- 
tograms trace the morphological distribution for the entire sam- 
ple. The dashed areas trace the histograms for the galaxies with 
velocity dispersions lower than 175 km s 



ter centers, where R200 is the radius within which the aver 
age cluster mass density is equal to 200 times the critica l den- 
sity. R ^on values were derived using equation (8) in Fin n et al 
(I2OO5I) and are listed in Table [1] These restrictions are gov- 
erned by the need to have the same spatial coverage in all 
structures given that a morphology-radius relation is observed 
in clusters (e.g. Whitmore & Gilmore 1991; Postman et al. 
2005). The fraction of galaxies from each morphological type 
have been corrected fo r magnitude and spatial incompleteness 



(iPoggianti et al.ll2006 ). Errors in the se fractions are calculated 



Gehrelsl (119861) . since Poissonian and bi- 



using the formulae in 
nomial statistics apply to our small samples. 



The fraction of early-type galaxies (E and SO) at redshifts 
z = 0.45, 0.55 and 0.75 are: 56!]^%, 61+^% and 75+^%, respec- 
tively. In other words, our calculations suggest that, in the red- 
sequence, the fraction of early-type galaxies decreases from 
z = 0.75 to z = 0.45. 

In order to check the robustness of this result, we also 
consider another technique of morphological classification. In 
particular the GIM2D decompositions presented in Simard et 
al. (2008, in preparation), where early-type galaxies are iden- 
tified by their bulge fraction (B/r> 0.35) and the ACS im- 
age smoothness within two half-light radii (52 < 0.075). This 



method leads to 33 



-I4c 



62!«'^ 



and 



80+^c 



of early-type 



galaxies at z = 0.45,0.55 and 0.75, respectively, in agreement 
with our previous finding. As mentioned earlier, only 3 clus- 
ters have ACS data in the 0.45 bin, and only 22 galaxies within 
O.6R200, as these clusters did not constitute our primary tar- 
gets for HST. This most likely explains the large difference 
found at this redshift between the two classification techniques. 
Nevertheless, the decrease of early-types fractions in the red- 
sequence is confirmed. It is worth noting that even when strong 
emission line galaxies are kept in the sample, the fractions of 
early-type galaxies are very much the same, i.e., 62^^^% and 
72+5% z=0.55 and z=0.75. 

The mean redsh ift of the MORPHS sample is 0.45 
(IDressler et al. I ll999!) . hence allowing direct comparison with 
our lowest redshift group and increasing the statistical sig- 
nificance (10 clusters). After selection of the red galaxies 
(g-r>1.4) without emission lines, i.e, a selection comparable to 
ours, the fraction of MORPHS Eh-SOs is 57+9% (B. Poggianti, 
private communication), i.e, very close to the initial value com- 
ing from our visual classification. Additional evidence of a 
decrease in early-type fraction w ith time on the red-sequence 
comes from lHolden et al.l ( 2006 ) who measured a total fraction 
of E/SO galaxies in two clusters at z=0.83 (MS 1054.4-0321 
and CI 0152.7-1357) of 74% and 82%, for their selection in 
mass ( > lO''''' Mq ), best corresponding to our red-sequence 
criteria. Note that our selection extends to lower masses or ac- 
cordingly to fainter limits, i.e., 3 - 6 x 10'" Mp (derived from 
stellar M/L with the calibration bvlBell & de Jond (1200 ih and 
the re-normalization proposed by ' de Jong & Belli "! 2007 ) for el- 
liptical galaxies). However, in this fainter regime, Holden et al. 
take into account all cluster galaxies, rather than only those 
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Fig. 14. Variation of the fraction of early-type galaxies in the 
red-sequence with redshift. Black and red symbols show the 
fraction derived from the EDisCS data. Black triangles show 
the results for the morphological classification based on visual 
inspection. Red squares stand for the G1M2D bulge and disk 
decomposition method, where early-type galaxies are selected 
upon their bulge-to-total light fractions and image smoothness. 
The green circle represents MORPHS and the blue star indi- 
cates the results of Holden et al. (2006), for their two X-ray 
selected clusters. 



on the red-sequence, making impossible a direct comparison. 
Anyhow, their figure 2 shows that their early type fraction can 
not vary much along the red-sequence. Figure [14] summarizes 
our results. 

The decrease, on the red-sequence, of the early-type frac- 
tion represe nts a remarkable synthesis between the results of 
Desai et al ] (&007) who detect no evolution in the cluster to- 
tal fraction of E/ SO galaxies b etween redshift 0.4 and 1 .4, and 
those of iDe Lucia et al. I (l2007l) who measure a decrease of the 
red-sequence faint-to-luminous galaxy ratio of ~40 % between 
redshift 0.45 and 0.75. I ndeed, if the red-seque nce keeps build- 
ing up, as suggested bv iDe Lucia et alj ( l2007h . it must find its 
supply in the cluster blue galaxies. These EDisCS blue galax- 
ies are composed of only 20+7% early-types and of 80% spi- 
ral galaxies. This distribution in types is constant with red- 
shift. Hence, provided that galaxy colours redden before any 
morphological transformation occurs, an increase in the num- 
ber of spirals on the red-sequence is naturally expected when 
new galaxies reach it. By restricting our analysis to the EDisCS 
spectroscopic sample, we do not reach as faint magnitudes as 



De Lucia et al. 



200 %. As a consequence, our 20% variation 
in the fraction of E/SO galaxies on the red-sequence between 
redshifts 0.45 and 0.75 most certainly constitutes a lower limit. 

It is interesting to estimate the subsequent evolution of 
the red-sequence population, i.e., from z - 0.45 to z = 0. 



Dressier et alJ (Il997h and lFasano et al. I (l2000h found that, while 
the elliptical fraction does not evolve in the redshift range 0.5 
to 0.0, the fraction of lenticular galaxies doubles. This increase 
in the number of SOs is accompanied by a drop in the frac- 
tion of (Sp-i-Ir r) type galaxies . Combining these res ults together 



with those of ^ Postman et al.l (12005 ) and EDisCS, iDesai et al 



12007) conclude that z ~ 0.4 must constitute a special epoch 
after which the total fraction of SO-type galaxies in clusters 
begins to increase. Meanwhile. iDe Lucia et al.l (120071) measure 
only ~ 10 % variation of the red-sequence faint-to-luminous 
galaxy ratio between redshift 0.45 and 0. As a consequence, at 
the time the red sequence reduces its growth (z < 0.45), one 
can expect a rise in its E/SO fraction as well. 

10. Discussion 

We have shown that the rate at which red-sequence galaxies 
evolve depends on their mass range. In particular, less mas- 
sive galaxies show evidence for a more extended star formation 
history. Two possibilities are considered: i) the sequence is in 
place at redshift z ~ 0.8, but a low le vel of star for mation con- 
tinues in galaxies on the faint end (e.g. lChiosi & Carraro 20021) : 
or ii) the red-sequence is continuously built-up by new incom- 
ers, preferentially selected among faint systems. We examine 
bot h hypotheses. 



Gebhardt et al.l (l2003h explore, in the context of the first 



model, the maximum amount of low-level star formation that 
the galaxies can experience without leaving the red-sequence. 
They found that, to keep the colours of the red-sequence galax- 
ies almost constant since z = 1, in agreement with observations, 
only 7%0 of the mass could have been formed after the initial 
burst (assumed to be at z =1.5) using exponential star forma- 
tion histories with e-folded time of 5 Gyr. We did a similar 
experiment, and calculated the maximum amount of low-level 
star formation that the galaxies could experience from z - 0.75 
to z = 0.45 in order to reproduce the evolution of the index-cr 
relations with redshift measured in Sect. [8] The exact percent- 
age depends on the chosen metallicity but we also find that this 
fraction has to be always less than ~10% (note that our redshift 
baseline is smaller). However, recent results of the evolution of 
the UV colou rs in early-type galaxies suggest that this percent- 
age is higher: iKavirail ( l2008h found that, at intermediate lumi- 
nosity (My < -20.5), early-type galaxies in the red-sequence 
have formed 30-60% of their mass since z = 1 . However, this 
study does not sample exclusively galaxies in dense clusters, 
but in a wide ra n ge of d i f Ferent environm e nts. 



Faber et alj (l2007l) : IPoggianti et"al] (l2006l) proposed 



mixed scenario where red-sequence ga laxies form through by 



two different channels. IPoggianti et al.l (| 2006 ) discussed a see 



nario where the red-sequence is composed of primordial galax- 
ies formed at Zf >2.5 on the one hand and of galaxies that 
quenched their star formation due to the dense environment of 
clusters on the other hand. They calculate that, of the 80% of 
passive galaxies at z-Q, 20% are primordia l and 60 % have been 
quenched. With regards to this later aspect, Harker et al. ( 20061) 



the percentage varies between 4 and 10% depending on the chosen 
metallicities 
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could reproduce the evolution of U — B and the H6, as well as 
the galaxy number density, by assuming models of constant star 
formation histories truncated at evenly spaced intervals of 250 
Myr up to now. 

Along this line, we calculate the number of new quenched 
galaxies that have to enter the red sequence to reproduce the re- 
sults obtained in Sec. 15. 21 We start with a population of galaxies 
with a formation redshift of z/ - 1 .4 and solar metallicity, both 
derived for the low-mass galaxy bin at z = 0.75. Then we con- 
sider galaxies that started forming stars at z = 2 at a constant 
rate of IMo/yr and start quenching them at regular intervals of 
250 Myr Once a galaxy is red enough to pass our red-sequence 
selection criteria at z = 0.45 and 0.55, we add its normalized 
spectrum to the that of the initial population. If we assume that 
all the new galaxies have velocity dispersions lower than 175 
km s"\ we need ~40% of new arrivals in the interval z=0.75 to 
z=0.45 in order to reproduce the constant luminosity-weighted 
age of the low-cr galaxy bins with redshift. Th is is in surpris 



ingly exact correspondence with the result of IPe Lucia et al 



(l2007h . who measure a decrease of the red-sequence faint-to- 
luminous galax y ratio of ~40 % betw een redshifts 0.45 and 
0.75. Note that IPe Lucia et al.l (|2007|) deal with photometric 



data and therefore their dataset reaches fainter magnitudes than 
our spectroscopic sample. The exercise presented here is of ex- 
tremely simple nature and is not intended to represent the true 
star formation histories of the galaxies. Apart form quenching, 
disk galaxies can merge and have bursts of star formation be- 
fore becoming red. By merging with other galaxies, our objects 
can also move from the low-cr bin to the high-cr bin, which is 
something not considered here. However, simple models simi- 
lar to this one have proved to be very successful in reproducing 
the evolution of the colours and Hd index and the evolution of 
the luminosity function (e.g.. Barker et al. 2007). 

Linking these constraints to our morphological results, it 
appears that the mechanism that quenches star formation does 
not necessarily produce a morphological transformation at the 
same time, but certainly, it favors it. Most likely, the mecha- 
nism that quenches star formation also produces the morpho- 
logical transformation, but on a different time-scale. Indeed, the 
fraction of spiral galaxies evolves by 20% and not by the 40% 
calculated in the spectroscopic analysis. Splitting between stel- 
lar population and morphologies was also raised by MORPHS 
dPressler et al.ll 19991 IPoggianti et alJll999l) . 



11. Summary 

We addressed the questions of the epoch of formation of the 
reddest galaxies in clusters, the extent of their period of star 
formation and the link between morphological and stellar pop- 
ulation evolution time scales. 

Our analyses are primarily based on 215 red-sequence 
galaxies, selected from the EDisCS spectroscopic database, 
that we divided into three redshift bins, z = 0.75, 0.55, and 
0.45. We considered their mass range, via the proxy of their 
velocity dispersions; their stellar population properties derived 
from absorption line features; and their morphologies, thanks 
to HST/ACS imaging. We have been able to trace the evolu- 
tion of the red-sequence from z = 0.75 to z = 0.45 in a ho- 



mogeneous dataset, the largest to date, extracted from a single 
survey, hence avoiding a mix of different systematic errors. 

• Before discarding red-sequence galaxies with strong 
emission lines from our absorption line analysis, we investi- 
gated the nature of their ionizing sources. Most of the EDisCS 
red-sequence galaxies with emission lines seem to be forming 
stars. The proportion of dusty star forming galaxies amon g our 
total s ample is larger than the fraction reported by iVan et al. 
(I2OO6I) at z = in the SDSS field red population. Whether this 



difference is due to evolution or to environment is worth clari- 
fying in the future. 

• We measured 12 Lick/IDS indices, carefully visually in- 
specting the spectra in order to eliminate those possibly af- 
fected by sky subtraction residuals or any other systematic ef- 
fect. To compare with our local sample of Coma galaxies, we 
derived new aperture corrections. We used state-of-the art stel- 
lar population models to derive ages, metallicities, and chemi- 
cal abundance ratios. These models make predictions over the 
whole spectral energy distribution, allowing us to analyze our 
data at a resolution of 325 km/s, avoiding any spurious correc- 
tion for the galaxy velocity dispersion. 

• After selecting on secure redshifts and signal-to-noise ra- 
tios, we stacked the galaxy spectra in redshift bins. In each 
bin, we also distinguished galaxies with respect to their veloc- 
ity dispersions, dividing at cr = 175 km s"'. We derived the 
age, metallicity (Z), and ff-element abundance ([E/Fe]) of each 
redshift and velocity dispersion group. 

- Massive galaxies ( cr > 175 km s"' ) show a variation in age 
corresponding to the expected cosmological variation be- 
tween redshifts, a mean solar metallicity and an overabun- 
dance [E/Fe] with respect to the solar values. Therefore, 
they are well represented by a scenario of formation at 
high redshift, followed by passive evolution. Conversely, 
the properties of less massive galaxies (cr < 175 kms~^) 
require longer star formation, at low level. Indeed, their lu- 
minosity weighted ages is found constant with time in our 
redshift range. An immediate consequence is that the age 
difference between low- and high-cr galaxies increases with 
time. This need to be taken into account in studies of the 
evolution of the colour-magnitude diagrams. 

- Values of [E/Fe] and [Z/H] are constant with time, inde- 
pendent of cr. This implies that the bulk of the galaxy stel- 
lar population is formed on a time-scale and with efficiency 
fixed by the galaxy mass. Possible subsequent episodes of 
star formation, which change the age of the less-massive 
galaxies, can not account for a large fraction of the galaxy 
mass. 

• We confirm that the evolution of the zero-points of 
the index-cr relationships with redshift is compatible with 
a scenario in which galaxies formed all their stars at high 
redshift and evolved passively since then. However, we 
demonstrate that it is also compatible with more complex star 
formation histories. In particular, galaxies can progressively 
enter the red-sequence as their star formation is quenched. In 
this scenario, galaxies with lower cr enter the red-sequence at 
lower redshift than more massive galaxies. 
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The two more important results of this work are: 

• The morphological analysis, based on visual or auto- 
mated classification, indicates that red sequences are composed 
of types covering the entire Hubble sequence, irrespective of 
the redshift, with, however, a very clear plurality of Es and SOs. 
The fraction of red-sequence early type galaxies decreases be- 
tween z - 0.75 and z - 0.45 by about 20%. This means that 
spiral galaxies get first redder, stopping their star formation, 
before they turn into earlier morphological types. This evolu- 
tion is however lower than the rate at which new galaxies enter 
the red-sequence, meaning that quenching and morphological 
transformation operate on different time-scales. 

• We found that the evolution of the line-strength indices 
with redshift can be reproduced if 40% of the galaxies with cr 
lower than 175 kms"' entered the red-sequence between z - 
0.75 and z - 0.45, in agreement with the fraction derived in 
studies of the luminosity function. 

In summary, a number of works dedicated to the evolution 
of the luminosity function of red galaxies in clusters (Rudnick 
et al. 2008, in preparation; De Lucia et al. 2004; De Lucia et 
al. 2007) and in the field (Faber et al. 2007; Brown et al. 2007) 
report an increase of the galaxy number density since z ~1. 
Conversely, colour-magnitude diagrams and index-cr relations 
of cluster red-sequence galaxies at intermediate- and high- 
redshift usually describe them as in place ver y early, passively 
evolving, independently of their mass (e.g., Ellis et alj 19971: 



Dressier et al 



1998 



2003 



1997 



Kelson et al. 
IWuvts et^lT 



van Dokkum et al 



19981 ; IStanford etail 
.2000 ; Zieg ler et al.1 12OOII; Blakeslee et al. 
2004 lTranetal.ll2007l) . The present work 



reconciles the luminosity functions and stellar content out- 
comes. 
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Appendix A: Transformation to the Lick system 

Although the new generation of stellar population models is 
starting to make predictions, not only of individual features, 
but of the whole spectral energy distributiori (e.g. IVazdekisI 
19991 IWMsMs.SLalJIaOOSUBruzual & Charlot"2003^, some of 
the most popular ones (e.g.|Worthey 1994; Thomas et al. 2003) 
still base their predictions on the original fitting functions re- 
leased by the Lick group dGorgas et al. 1993 : Worthey et al 



1994|). In order to compare with these models and with the bulk 
of the data in the literature, which is usually transformed onto 
this system, one has to follow two steps. The first step is to 
degrade the resolution of the spect ra into the wavelength vari- 
able resolution of the Lick system ( Worthey & Ottavianill 19971 
see) and then apply a correction due to the velocity dispersion 
of the galaxies. The second step is to correct small system- 
atic differences in the indices due to differences in the shape 
of the spectra of the Lick/IDS stars and a fuUy-flux calibrated 
spectrum. To do that, one usually observes stars in common 
with the Lick library with the same instrumental configuration 
as the objects one wants to analyze and compare the indices 
after degrading their spectra to the Lick resolution. This can- 
not be done, though, in the growing number of work analyz- 
ing Lick indices in galaxies at high redshift. For these cases 
when stars in common with Lick are not present, it is still pos- 
sible to correct for the instrumental response of the Lick/IDS 
stars by applying an offset derived by comparing flux cali- 
brated spectra with Lick/IDS spectra previously degraded to 
the Lick/I DS resolution. This experim ent has been attempted 



before bv lWorthev & Ottavianil ( 1 19971) andlBruzual & Chariot! 
(l2003h using the Jones (1999) and STELIBlLe Borgne et al. 



(120031) 

2003 ') ] 



2003) libraries respe ctively. Here we present new offsets using 
the MILES library (iSanchez-Blazquez et alJl2006c ). The ad- 
vantage of this calibration is that the MILES library contains 
a large number of stars (242) in common with the Lick/IDs li- 
brary. Furthermore, the flux calibration is very accurate (better 
than 1%) 

To obtain the offsets we broadened the MILES spectra to 
the wavelength-d ependent resolution of the L ick/IDS spectra 



using figure 7 in Worthey & Ottavianil (1997) a nd measured 



Index 


offset 


Linear relation 




/^A/TTT T 

(^IVllLno-LlCK ) 


(MILES =a+ b x Lick) 


CNi 


0.007 


0.0082(0.0017) 


+ 0.9112(0.0124) 


CN2 


0.006 


0.0101(0.0021) 


+ 0.9040(0.0136) 


Ca4227 


-0.081 


0.0627(0.0256) 


+ 0.9046(0.0119) 


G4300 


-0.032 


0.5795(0.0663) 


+ 0.8480(0.0126) 


Fe4383 


-0.587 


-0.1660(0.0798) 


+ 0.9321(0.0154) 


Ca4455 


-0.391 


-0.0752(0.0340) 


+ 0.7849(0.0204) 


Fe4531 


-0.109 


0.3201(0.0544) 


+ 0.8541(0.0164) 


Fe4668 


0.050 


0.5664(0.0649) 


+ 0.8873(0.0110) 


HyS 


0.068 


0.1535(0.0283) 


+ 0.9755(0.0108) 


Fe5015 


-0.629 


0.1940(0.0543) 


+ 0.8573(0.0088) 


Mgi 


-0.005 


-0.0024(0.0001) 


+ 0.9263(0.0082) 


Mg2 


-0.020 


-0.0043(0.0012) 


+ 0.9186(0.0050) 


Mgb 


-0.208 


0.2201(0.0351) 


+ 0.8934(0.0089) 


Fe5270 


-0.151 


0.1311(0.0393) 


+ 0.9008(0.0147) 


Fe5335 


-0.060 


0.1059(0.0361) 


+ 0.9268(0.0150) 


Fe5406 


0.045 


0.2190(0.0242) 


+ 0.8543(0.0151) 


Fe5709 


-0.029 


0.0135(0.0227) 


+ 0.9162(0.0242) 


Fe5782 


-0.029 


0.0157(0.0208) 


+ 0.9361(0.0267) 


Na5895 


-0.148 


-0.1842(0.0352) 


+ 1.0038(0.0107) 


TiO, 


-0.014 


-0.0066(0.0006) 


+ 0.9128(0.0065) 


TiOz 


-0.010 


0.0015(0.0007) 


+ 0.9119(0.0045) 


uSa 


0.083 


0.0286(0.0512) 


+ 0.9354(0.0115) 




-0.073 


-0.2029(0.0565) 


+ 0.9607(0.0082) 




-0.013 


0.0069(0.0329) 


+ 0.9585(0.0147) 




-0.051 


-0.0750(0.0258) 


+ 0.9763(0.0091) 



Table A.l. Mean index offsets between he stars in common be- 
tween MILES and the Lick/IDS. Some times the differences 
can be better described with a linear relation than a simple off- 
set. Last column display the best linear relation describing the 
relation of the indices measured in MILES and Lick/IDS spec- 
tra. 



than with a single offset. Therefore, we have also derived a lin- 
ear transformation between both libraries obtained with a sim- 
ple linear regression eliminating galaxies deviating more than 
3 times the rms from the relation. 

In principle, the offsets between the indices measured in the 
Lick and other libraries are due to differences in the continuum 



shape. However, as was noted before by IWorthev & Ottaviani 
(1997), some of the indices with larger differences have very 
narrow definition bands and show a very small dependence 
with the flux shape. Therefore, the reason for the differences 
is not clear. 

The mean offsets and linear relations obtained in the com- 
parison between the MILES and Lick/IDS libraries are listed in 
TablelAH 



the L ick/IDS indices with the definition given in iTrager et al 
(fl998,) . In some cases, differences between the indices mea- 
sured in both datasets can be described with a linear fit, better 



Appendix B: IVIean gradients and aperture 
correction 

Galaxies have radial gradients in the line-strength indices and, 
therefore, when comparing galaxies at different redshifts , the 
indices need to be aperture corrected. |j0rgensenl (119970 and 
J0rgensen et alj ( 120051) pr esent a series of correc t ions based on 



the rad ial gradients from IVazdekis et alj ( 11997b . ICardiel et al . 
(Il998l) . and lDavidge & Clarkl(ll994 . However, fliose works do 
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CN1 CN2 Ca4227 G4300 Fe4383 Ca4455 




5 10 

Lick 



Fig. A.l. Comparison between the Lick/IDS indices measured in the stars in common between the MILES and Lick/IDS Ubraries. 
The soUd hne represents the 1 : 1 relation. The dashed line shows the calculated offset, while the red dashed line indicates the best 
linear fit to the data. 



not include many of the Lick indices, particularly the higher 
order Balmer lines. 

We have derived new aperture corrections for all the 
Lick/IDS indices + the hi gher order Balmer line indices 
( Worthey & Ottaviaiiil 1997 ). taking advantage of the large 
sample of ne arby galaxies presented in Sanchez-Blazquez et al 



(l2006ah and ISanchez-Blazquez et al. I (l2007h . To calculate the 



a parameter, we assume that the light in our galaxies follows 
a de Vaucouleur profile. Using the gradients for each individ- 
ual galaxy, we obtain the final index inside an aperture as the 
luminosity-weighted index. We then calculate how the final in- 
dices depends on aperture size. The variation of the luminosity- 
weighted indices with different aperture sizes can be described 
as power law. Therefore, aperture corrections (those measured 



in A) can be written as (see |j0rgensen et alJll995l : |j0rgensen 



1997D: 



log{index)corr — \og(index)ap + a log (B.l) 

^nor 

and for molecular indices and higher-order Balmer line^ 
indeXcorr — indeXgp + a log — ^, (B-2) 

^ nor 

where rgp is the equivalent circular aperture, obtained as 2rap - 
1.025 X 2{xyln)^l^; and x and y are the width and the length of 
the rectangular aperture. 



' The correction for the higher-order Balmer lines is additive as 
these indices can have both positive and negative values. 
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Index 


O! 


Type 


cr 


0.058 ± 0.047 


1 


D4000 


0.027 ±0.018 


1 


H<5a 


-0.721 ±0.508 


2 


mp 


-0.206 ± 0.249 


2 


CNi 


0.058 ± 0.023 


2 


CN2 


0.059 ± 0.030 


2 


Ca4227 


0.047 ± 0.075 


1 


G4300 


0.027 ± 0.022 


1 


HyA 


-0.850 ± 0.422 


2 


H7f 


-0.436 ± 0.321 


2 


Fe4383 


0.061 ± 0.033 


1 


Ca4455 


0.097 ± 0.075 


1 


Fe4531 


0.040 ±0.021 


1 


Fe4668 


0.119 ±0.049 


1 


Fe5015 


0.05 ± 0.02 


1 


W 


-0.05 ±0.19 


1 


Mgb 


0.048 ±0.051 


1 



Table B.l. a coefficient and corresponding error obtained using 
the samples of SB06 and Sanchez-Blazquez et al. 2007. 



The a coefficients are listed in Table IbTI and can be used 
by any other study. In our case, we were more interested in 
selecting a fraction of the total light of the galaxy instead of a 
fixed aperture and therefore y is not fixed but is chosen to be 
a FWHM of the spatial profile. However, we have fixed the y 
for all the spectra to be the physical aperture equivalent to 1 
arcsec (the width of the slit) at z = 0.75. The correction for the 
galaxies at z > 0.4 was very small and we did not apply any. 
The aperture correction was applied only to the galaxies in the 
Coma cluster 

Elliptical galaxies show a great variety of gradients that 
do not seem to be correlated with any other properties of the 
galaxies, like mass or colour For this reason, we used the mean 
gradient to compute the aperture corrections. We calculated an 
error due to the dispersion among the mean of all the mea- 
sured gradients. To do this, we performed a series of Monte 
Carlo simulations where each gradient was perturbed follow- 
ing a Gaussian distribution with cr equal to the typical RMS 
dispersion in the gradients. For each new gradient, the coeffi- 
cient a and the final correction for the galaxies in the Coma 
cluster were recomputed. The errors in the aperture correction 
calculated this way are indicated in TablelH 

Appendix C: Relation of tlie Lick/IDS indices with 
velocity dispersion 

Figure [TO] shows the relation between all the measured line- 
strength indices and the velocity dispersion for the 4 considered 
redshift bins. 
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Fig. C.l. Relation of the indices with the velocity dispersion for the galaxies of the sample grouped in three different redshift 
bins. The meaning of the lines is the same as in Fig. [10] From top to bottom: z - 0.75, z - 0.55, z - 0.45, and Coma galaxies are 
represented. 
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Group Cluster Name 



ACS N,.,, 



0.45 



0.75 



a ± err 
(kms-') 



C11018.2-1211 

cll037.9-1243a 

C11059.2-1253 

C11138.2-1133 

clll38.2-1133a 

C11202.7-1224 

C11301.7-1139 

cll301.7-1139a 

C11420.3-1236 



0.4734 
0.4252 
0.4564 
0.4796 
0.4548 
0.4240 
0.4828 
0.3969 
0.4962 



17 
18 
26 
17 
8 

18 
19 
13 
18 



9 
9 

18 
8 
3 
7 

10 
6 

13 



CQ7+46 

510+'^ 

•732+72 

-76 

542+5] 
518-104 

218+« 



C11040.7-1155 

C11054.4-1146 

C11054.7-1245 

clll03.7-1245b 

C11216.8-1201 

C11354.2-1230 



0.7043 
0.6972 
0.7498 
0.7031 
0.7943 
0.7620 



14 
36 
20 

3 
38 

7 



8 418!: 
16 

6 

2 
23 

5 



46 

589!?g 
504!^^3 

252!« 
1018^^3 

648+!?^ 



R200 
Mpc 



0.93 
1.06 
0.99 
1.40 
1.05 
1.02 
1.31 
0.78 
0.41 



0.55 C11037.9-1243 


0.5783 


X 


7 


5 


319!M 


0.57 


clll03.7-1245a 


0.6261 


X 


9 


6 


336!1« 


0.59 


C11119.3-1129 


0.5500 




17 


9 


166!?^ 


0.30 


C11227.9-1138 


0.6357 


X 


14 


8 


574+'2 


1.00 


cll227.9-1138a 


0.5826 


X 


4 


1 


341+**^ 


0.61 


C11232.5-1250 


0.5414 


X 


41 


20 


1080+^^' 


1.99 


C11353.0-1137 


0.5882 




10 


7 


666!|3« 


1.19 


cll354.2-1230a 


0.5952 


X 


8 


5 




0.77 


C11411. 1-1148 


0.5195 




18 


12 




1.32 



0.70 
1.06 
0.70 
0.42 
1.61 
1.05 



Table 1. Col. (1) Redshift in terval in which the galaxy structure falls (see text for details). Col. (2): Cluster name. Following 
Milvang- Jensen et al.l (120081) . we label the secondary structures with letters "a" and "b". Col. (3): Cluster redshift. Col. (4): 
marked with x if ACS images are available. Col. (5): Number of secure spectroscopically confirmed red-sequence galaxies. Col. 
(6) Number of red-sequence galaxies with no or negligible e mission, full-filli ng our signal-to-noise ratio criterion. Col. (7) 
Velocity dispersion of the cluster. Col. (8) R200, calculated as in iFinn et alj (l2005b . 



Table 3. Line-strength indices in our sample of red-galaxies with S/N(A)>10 and no- or weak-emission lines. The number 
between brackets after the indices values indicates the presence (1) or not (0) of sky-subtraction residuals or telluric absorptions 
inside the definition passband of the index. Only those indices with labels (0) are used in our analysis. Type: Indicate the presence 
(2) or not (1) of weak ([OII]<7A) emission lines. Last column lists the measured velocity dispersion and its error This is only a 
portion of the Table, shown for guidance regarding its format and content. The full table is electronically available at CDS and at 
^http://www. ucm.es/inf o/Astrof/psb/ediscs. html. 





D4000 


H<5a 


H(5f 


CN2 


Ca4227 


G4300 




1018471-1210513 


1.83±0.01 (1) 


1.28±0.41 (0) 


2.11±0.27 (0) 


0.027±0.015 (0) 


1.27= 


tO.23 (0) 


2.71 ±0.43 (0) 


-1.53±0.47 (0) 


1018464-1211205 


2.22±0.02 (1) 


-0.18±0.69 (0) 


0.84±0.47 (0) 


0.038±0.024 (0) 


0.87= 


tO.38 (0) 


4.86±0.65 (0) 


-4.75±0.81 (0) 


1018467-1211527 


2.13±0.01 (0) 


-0.97±0.29 (0) 


0.63±0.19(0) 


0.108±0.010(0) 


0.48= 


t0.16(0) 


5.10±0.26 (0) 


-5.55±0.33 (0) 


1018401-1214013 


1.98±0.02 (0) 


-0.57±0.65 (0) 


0.40±0.45 (0) 


0.072±0.022 (0) 


0.53= 


tO.37 (0) 


7.13±0.57 (0) 


-4.71 ±0.76 (0) 






Fe4383 


Ca4455 


Fe4531 




type 


cr±err 




1018471-1210513 


0.52±0.28 (0) 


3.71±0.68 (0) 


0.67±0.38 (0) 


3.28+0.58 (0) 




2 


245.9±15.8 




1018464-1211205 


-2.01±0.52 (0) 


4.10±1.09 (0) 


0.91 ±0.60 (0) 


3.47±0.93 (0) 




1 


279.9±27.6 




1018467-1211527 


-1.73±0.21 (0) 


3.27±0.44 (0) 


0.98±0.24 (0) 


3.03±0.38 (0) 




1 


211.1±14.8 




1018401-1214013 


-0.29±0.45 (0) 


1.77±1.06 (0) 


0.71±0.57 (0) 


1.37±0.91 (0) 




1 


129.2±21.0 
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Index 


Coma 


■7 — n AS 

Z — U.^J 


■7 — n 
z — y}.jj 


■7 — 

Z — yJ. Jj 


correcuon lype 




1.031 


0.063 


0.016 


0.00 


additive 


CNi 


-0.083 


-0.005 


-0.001 


0.00 


additive 


Ca4227 


0.856 


0.990 


0.997 


1.0 


multiplicative 


G4300 


0.914 


0.995 


0.998 


1.0 


multiplicative 


Hy,- 


0.623 


0.038 


0.009 


0.0 


additiva 


Fe4383 


0.817 


0.987 


0.996 


1.0 


multiplicative 


Ca4455 


0.726 


0.980 


0.995 


1.0 


multiplicative 


Fe4531 


0.876 


0.992 


0.998 


1.0 


multiplicative 


C4668 


0.676 


0.976 


0.994 


1.0 


multiplicative 


H/3 


1.179 


1.010 


1.000 


1.0 


multiplicative 



Table 4. Aperture correction for Lick/IDS indices measured in our sample. The correction for the molecular indices and for the 

higher-order Balmer lines is additive while the correction for the atomic indices is multiplicative. The aperture correction was 
only applied to the Coma galaxies, as it is very small for the other redshifts. The type of correction is indicated in the last column. 



Redshift group 






D4000 


USa 


H(5f 


CN2 


Ca4227 


G4300 


HyA 


Hyf 


Fe4383 


cr 








A 


A 


A 


mag 


A 


A 


A 


A 


A 


kms"' 


0.45 


cr > 


175 km/s 


2.044 


-0.433 


0.894 


0.067 


0.932 


4.589 


-4.114 


-0.704 


3.662 


233.1 








0.004 


0.140 


0.094 


0.005 


0.076 


0.133 


0.161 


0.098 


0.220 


7.2 


0.45 


cr < 


175km/s 


2.045 


0.103 


1.074 


0.036 


0.841 


4.053 


-3.145 


-0.177 


3.380 


122.3 








0.004 


0.139 


0.094 


0.005 


0.077 


0.136 


0.158 


0.096 


0.222 


12.0 


0.55 


cr > 


175 km/s 


2.051 


-0.075 


1.143 


0.065 


0.847 


4.600 


-3.779 


-0.437 


3.534 


218.1 








0.002 


0.087 


0.058 


0.003 


0.048 


0.083 


0.100 


0.061 


0.137 


9.0 


0.55 


cr < 


175 km/s 


1.995 


0.045 


0.960 


0.028 


0.757 


4.341 


-3.138 


-0.114 


3.138 


134.9 








0.004 


0.140 


0.095 


0.005 


0.078 


0.135 


0.159 


0.096 


0.224 


10.1 


0.75 


cr > 


175 km/s 


2.008 


0.699 


1.400 


0.057 


0.824 


4.381 


-3.061 


-0.234 


2.877 


233.2 








0.003 


0.104 


0.070 


0.004 


0.059 


0.102 


0.120 


0.073 


0.171 


6.2 


0.75 


cr < 


175 km/s 


1.987 


0.453 


1.163 


0.034 


0.826 


4.244 


-2.963 


-0.016 


2.796 


154.2 








0.003 


0.130 


0.089 


0.004 


0.073 


0.128 


0.149 


0.091 


0.213 


8.1 



Table 5. Line-strength indices and velocity dispersion measured in the stacked spectra. The second row of each line Usts the 
errors. Last column show the velocity dispersion measured in the different stacked spectra. 







cr > 175 km s ' 






cr < 175 km s ' 






Age 


[Z/H] 


[E/Fe] 


Age 


[Z/H] 


[E/Fe] 




Gyr 






Gyr 






Coma 


7.90 ± 1.81 


-0.03 ± 0.06 


0.40 + 0.03 


5.80 + 0.36 


-0.11 + 0.03 


0.31+0.04 


z = 0.45 


5.60+1.05 


-0.01 ± 0.07 


0.35 + 0.09 


3.20 + 0.64 


-0.05+0.13 


0.16 + 0.12 


z = 0.55 


3.40 + 0.52 


0.11 +0.07 


0.35 + 0.06 


3.40 + 0.81 


-0.13 + 0.09 


0.11+0.13 


z = 0.75 


2.70 + 0.38 


0.18 + 0.11 


0.46 + 0.06 


3.10 + 0.83 


-0.08 + 0.14 


0.30 + 0.21 


z = 0.45 


5.70 + 2.42 


0.02 + 0.12 


0.45 + 0.10 


2.80 + 0.46 


0.01+0.11 


0.23 + 0.14 


z = 0.55 


3.10 ±0.62 


0.11+0.12 


0.30 + 0.10 


2.90 + 0.42 


-0.04 + 0.10 


0.25+0.11 


z = 0.75 


4.60 + 2.32 


-0.07 + 0.14 


0.55 + 0.10 


2.50 ± 0.34 


0.11+0.13 


0.24 + 0.15 



Table 6. First 4 rows: Ages (measured in Gyr), metallicities [Z/H] and [E/Fe] abundances for the stacked spectra, in different 
redshift bins, of the red-sequence galaxies, measured using V08 models and the indices Hy^ , Fe4383 and CN2. The spectra have 
been stacked separating galaxies with cr lower than 175 km s"' (right side) and larger than 175 km s"^ (left side). The parameters 
have been obtained iteratively as described in the text. Last three rows: Ages, metallicities and [E/Fe] derived for the stacked 
spectra, at different redshifts, of morphologically classified E and SO galaxies. 
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linear fit 




Zf = 1.4 




Zf - 2 




Zf - 3 




Index 


reusnitt (S'ZUU Km/s 


(T 


^Zuv Km/s 


t 


IS'ZUU Km/s 


t 


t^PZUU Km/s 


t 


Fe4383 


0.02 4.70 


0.41 


4.09 


1.5 


4.33 


0.9 


4.38 


0.8 




0.45 3.58 


0.82 


3.39 


0.2 


3.57 


0.0 


3.65 


0.1 




0.55 3.73 


0.38 


3.06 


1.7 


3.40 


0.8 


3.47 


0.8 




0.75 2.77 


1.24 


2.56 


0.2 


3.19 


0.3 


3.34 


0.4 


USA 


0.02 -1.40 


0.79 


-1.37 


0.0 


-1.85 


0.6 


-1.93 


0.7 




0.45 -0.23 


1.12 


-0.08 


0.1 


-0.47 


0.2 


-0.64 


0.4 




0.55 0.09 


1.09 


0.58 


0.4 


-0.02 


0.1 


-0.26 


0.3 




0.75 0.63 


1.10 


1.54 


0.8 


0.33 


0.3 


0.05 


0.5 


HyF 


0.02 -1.48 


0.44 


-0.86 


1.4 


-1.09 


0.9 


-1.13 


0.8 




0.45 -0.61 


0.66 


0.10 


1.1 


-0.36 


0.4 


-0.48 


0.2 




0.55 -0.48 


0.66 


0.31 


1.1 


-0.09 


0.6 


-0.21 


0.4 




0.75 0.04 


0.97 


0.96 


0.9 


0.16 


0.12 


-0.04 


0.0 



Table 7. Comparison between the indices predicted by the linear fit to the data at cr = 200 km s and the ones predicted - at the 
same velocity dispersion- in a passively evolving model with redshif t formations zi = lA,Zf = 2 and Zf = 3 . Column 3 shows the 
predicted value by the linear fit at cr = 200 km s~^ , column 4 the standard deviation among the relation; columns 5, 7 and 9 show 
the value predicted by the model with redshift formations of 1 .4, 2 and 3 respectively. Columns 6, 8 and 10 shows the f-parameter 
of comparing these two values (fitted and predicted) using a Student's t-test. A t value higher than 1.96 would indicate that the 
probability that the two values are different by chance is less than 5% (for samples larger than 40 and 6% for smaller samples). 
As can be seen, none of the values show significant differences from the model. 







linear fit 




Zf = 1.4 




Zf = 2 




Zf = 3 




Index 


redshift 


@200 km s"' 


cr 


@200 km s"' 


t 


@200kms-l 


t 


@200kms"' 


t 


D4000 


0.45 


2.03 


0.11 


1.94 


0.8 


1.98 


0.4 


2.01 


0.2 




0.55 


2.08 


0.07 


1.85 


3.2 


1.93 


2.1 


1.96 


1.7 




0.75 


1.98 


0.12 


1.75 


1.9 


1.89 


0.7 


1.92 


0.5 




0.02 


2.12 


0.12 


2.11 


0.1 


2.16 


0.3 


2.20 


0.7 


KSf 


0.45 


0.92 


0.52 


1.05 


0.2 


0.86 


0.1 


0.79 


0.2 




0.55 


0.99 


0.89 


1.36 


0.4 


1.07 


0.1 


0.97 


0.0 




0.75 


1.40 


0.71 


1.80 


0.6 


1.24 


0.2 


1.10 


0.4 




0.02 


0.70 


0.33 


0.50 


0.6 


0.38 


0.9 


0.28 


1.3 


Ca4227 


0.45 


0.82 


0.31 


1.06 


0.8 


1.12 


0.9 


1.15 


1.1 




0.55 


0.89 


0.27 


0.97 


0.3 


1.05 


0.6 


1.08 


0.7 




0.75 


0.83 


0.32 


0.86 


0.1 


0.99 


0.5 


1.05 


0.7 




0.02 


1.17 


0.16 


1.27 


0.6 


1.33 


1.0 


1.39 


1.3 


G4300 


0.45 


4.58 


0.81 


3.87 


0.9 


4.03 


0.7 


4.12 


0.6 




0.55 


4.41 


0.53 


3.55 


1.6 


3.85 


1.0 


3.95 


0.8 




0.75 


4.06 


1.31 


3.05 


0.8 


3.66 


0.3 


3.81 


0.2 




0.02 


5.25 


0.51 


4.35 


1.7 


4.41 


1.6 


4.45 


1.5 


Ca4455 


0.45 


0.86 


0.28 


0.88 


0.1 


0.91 


0.2 


0.91 


0.2 




0.55 


0.96 


0.30 


0.83 


0.4 


0.88 


0.3 


0.90 


0.2 




0.75 


0.79 


0.61 


0.75 


0.1 


0.85 


0.1 


0.87 


0.1 




0.02 


1.15 


0.12 


0.97 


1.5 


1.00 


1.2 


1.02 


1.1 


Fe4531 


0.45 


2.77 


0.43 


2.76 


0.0 


2.82 


0.1 


2.83 


0.1 




0.55 


2.93 


1.74 


2.66 


0.1 


2.75 


0.1 


2.78 


0.1 




0.75 


2.51 


0.80 


2.53 


0.0 


2.69 


0.2 


2.74 


0.3 




0.02 


3.04 


0.23 


2.92 


0.5 


2.98 


0.3 


3.03 


0.0 




0.45 


-3.85 


1.35 


-2.92 


0.7 


-3.38 


0.3 


-3.61 


0.2 




0.55 


-3.83 


1.77 


-2.12 


0.9 


-2.89 


0.7 


-3.12 


0.4 




0.75 


-2.56 


1.87 


-0.86 


0.9 


-2.41 


0.1 


-2.78 


0.1 




0.02 


-5.56 


1.04 


-4.38 


1.1 


-4.68 


0.8 


-4.90 


0.6 



Table C.l. Comparison between the indices predicted by the linear fit to the data at cr = 200 km s and the ones predicted - at 
the same velocity dispersion - in a passively evolving model with formation redshifts of Zf = 1.4, Zf = 2 and Zf = 3. Column 3 
shows the value predicted by the linear fit at cr = 200 km s Column 4 shows the standard deviation of the relation; ColumnS 
shows the value predicted by the SSP model; and the last column shows the t parameter of the comparison of these two values. 
A t value higher than 1 .9 would indicate that the there is a 5% probability that the two values are different by chance. As can be 
seen, none of the values show significant differences from the model. 



